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ABSTRACT 

A theoretical  approach  Is  developed  and  a computational  procedure  adapt- 
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and  the  effect  of  the  •'’’flow  distortion  thickness"  are  studied  by  means  of  the 
"thin  body"  approximation. 

The  surface  integral  equation  which  relates  the  unknown  loading  to  the 
known  velocity  distribution  on  the  blades  is  solved  by  the  mode  approach  in 
conjunction  with  the  "Tift  operator"l''technique.  The  analysis  treats  both  design 
and  off-design  conditions  in  steady-state  and  unsteady  flows,  and  the  proper 
chordwise  modes  are  selected  for  each  condition. 

The  numerical  solution  yields  the  blade  loading  and  resulting  hydrodynamic 
forces  and  moments  and  blade  bending  moments,  and,  in  addition,  the  blade 
pressure  distributions  on  each  blade  face  due  to  both  loading  and  thickness 
effects,  thus  providing  information  necessary  for  the  prediction  of  cavitation 
inception.  Calculations  have  been  performed  for  a set  of  three  3“bladed 
propellers  of  different  EAR  operating  in  a screen-generated  wake,  for  compar- 
ison with  experimental  data. 
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ABSTRACT 


A theoretical  approach  is  developed  and  a computational  procedure 
adaptable  to  a high  speed  digital  computer  is  established  for  the  evalua- 
tion of  the  blade  pressure  distribution  of  a marine  propeller  due  to 
thickness  and  loading  effects.  The  dual  role  of  the  blade  thickness  is 
considered.  The  contribution  of  the  "nonplanar  thickness"  to  the  propeller 
loading  and  pressure  distribution  and  the  effect  of  the  "flow  distortion 
thickness"  are  studied  by  means  of  the  "thin  body"  approximation. 

The  surface  integral  equation  which  relates  the  unknown  loading  to 
the  known  velocity  distribution  on  the  blades  is  solved  by  the  mode  ap- 
proach in  conjunction  with  the  "lift  operator"  technique.  The  analysis 
treats  both  design  and  off-design  conditions  in  steady-state  and  unsteady 
flows,  and  the  proper  chordwise  modes  are  selected  for  each  condition. 

The  numerical  solution  yields  the  blade  loading  and  resulting 
hydrodynamic  forces  and  moments  and  blade  bending  moments,  and,  in  addition 
the  blade  pressure  distributions  on  each  blade  face  due  to  both  loading 
and  thickness  effects,  thus  providing  information  necessary  for  the  pre- 
diction of  cavitation  inception.  Calculations  have  been  performed  for  a 
set  of  three  3-bladed  propellers  of  different  EAR  operating  in  a screen- 
generated wake,  for  comparison  with  experimental  data. 
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NOMENCLATURE 


A (A) 

function  defined  in  Eqs.  (33)  and  (35) 

a 

u/J  (nondimens ional ) inversely  proportional  to  advance 
ratio  at  design  condition 

aod 

7f/dQd  (nond imens ional ) inversely  proportional  to  advance 
ratio  of  off-design  condition 

a 

designation  of  NACA-a  meanline 

CF,CD 

frictional  coefficients 

CP 

pressure  coefficient 

c 

expanded  chord  length,  ft. 

F 

x,y,z 

propel ler- induced  forces  in  x,y,z  direction  (see  Fig.  1) 

f 

c 

camberline  ordinates  from  face  pitch  line 

f 

T 

blade  thickness  distribution 

l(  ) 

modified  Bessel  function  of  first  kind,  of  order  m 

, (m)  ( ) 

defined  in  Eq.  (20) 

i 

index 

J 

design  advance  ratio 

Jod 

off-design  advance  ratio 

AJ 

J j-J 
od 

j 

i ndex 

K 

kernel  of  integral  equation 

modified  kernel  after  chordwise  integrations 

M ) 

m 

modified  Bessel  function  of  second  kind  of  order  m 

k 

variable  of  integration 

Ltq) <0,ea) 

loading  distribution  in  Ib/ft. 

\ 
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L(ql  (r) 

spanwise  loading  in  lb/ft. 

♦ 

(p) 

spanwise  loading  coefficients  of  the  chordwise  modes  in  lb/ft 
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integer  multiple 

blade  bending  moment  about  face  pi  ten  line 

source  strength  at  point  U,p,Gq) 

index  of  summation 

order  of  lift  operator 

number  of  blades 

blade  index 

rps 

order  of  chordwise  mode 

unit  normal  vector  on  nelicoidal  surface  at  loading  point 

unit  normal  vector  on  helicoidal  surface  at  control  point 

pressure,  lb/ft2 

P_  - P+,  pressure  jump,  lb/ft2 

geometric  pitch  at  each  radial  position,  ft. 

propel  1 er- induced  moments  about  x,y,z  axis  (see  Fig.  1) 

order  of  harmonic  of  inflow  £ield 

Descartes  distance 

Reynolds  number  based  on  cnora  lengtn 

radial  ordinate  of  control  point 

propel ler  radius,  ft. 

propeller  lifting  surface,  ft2 

chordwise  location  as  fraction  of  chord  length 

time,  sec. 

maximum  thickness  of  blade,  ft. 
free  stream  velocity,  ft/sec.  (design) 
free  stream  velocity,  ft/sec.  (off  design) 
variable  of  integration 
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Fourier  coefficients  of  the  known  downwash  velocity 
distribution 

normal  velocity  due  to  camber  effects 
normal  velocity  due  to  flow-incidence  angle 
normal  velocity  due  to  wake 

normal  velocity  due  to  nonplanar  thickness 

longitudinal  perturbation  velocity  (V^-U) 

tangential  perturbation  velocity 

measured  axial  velocity 

velocity  distribution  normal  to  propeller 

defined  in  Equation  (20) 

longitudinal  ordinate  of  control  point 

cylindrical  coordinates  of  control  point 

Cartesian  coordinate  system 

tan  ' (u/flr)  = tan  1 (l/ar) 

defined  in  Eq.  (12) 

chordwise  mode  shapes 

-fit  (see  Fig.  2) 

angular  ordinate  of  loading  point 

angular  position  of  loading  point  with  respect  to  olade 
reference  line,  in  moving  coordinate  system  (see  Fig.  2) 

subtended  angle  of  projected  blade  semichord,  radians 

2ir (n-1  )/N , n»l,2,...N 

geometric  pitch  angle  at  each  radial  position 

angular  chordwise  location  of  loading  point  (in 
trigonometric  transformation  Eq.  14) 

defined  in  Eq.  (21)  (see  Appendix  A) 
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A|n)  (y) 


defined  in  Appendix  A 
kinematic  v i scos i ty 

longitudinal  ordinate  of  loading  point 

cylindrical  coordinates  of  loading  point 

radial  ordinate  of  loading  point 

mass  density  of  fluid 

ratio  of  leading  edge  radius  to  chord 

angular  measure  of  skewness  from  blade  reference  line  (see  Fig.  2) 
o'"  - op  = difference  in  skewness  at  control  and  loading  point 
variable  of  integration 
velocity  potential 

generalized  lift  operator  (Eq.  (18)) 
angular  ordinate  of  control  point 

angular  position  of  control  point  with  respect  to  blaoe 
reference  line,  in  mov ing* coord inate  system  (see  Fig.  2) 

angular  chordwise  location  of  control  point  (in  trigono- 
metric transformation  Eq.  15) 

acceleration  potential 

magnitude  of  angular  velocity  of  propeller 
angular  frequency  of  loading 
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INTRODUCTION 

Davidson  Laboratory  has  been  engaged  in  a series  of  investigations 
concerned  with  the  adaptation  of  the  linearized  unsteady  lifting  surface 
theory  to  the  case  of  a marine  propeller  operating  in  a nonuniform  flow 
field  (see,  for  example,  References  l-3)>  the  objective  of  which  has  been 
the  prediction  of  blade  loading  distributions,  propeller-generated  forces 
and  moments  and  blade  bending  moments.  In  these  investigations,  it  has 
been  assumed  that  the  lifting  surfaces  (blades)  have  zero  thickness. 

The  linearized  wing  theory  considers  that  an  arbitrary  wing  cross- 
section  is  composed  of  a symmetrical  part  corresponding  to  its  thickness 
distribution  and  an  asymmetrical  part  of  zero  thickness  but  with  camber 
and  angle  of  attack.  The  asymmetrical  part  contributes  to  the  lift  of 
the  wing.  On  the  other  hand,  the  thickness  distribution  contributes 
nothing  to  the  lifting  properties  of  a wing  unless  it  is  nonplanar. 

The  marine  propeller  with  its  blade  describing  a helicoidal  surface 
is  one  of  the  few  nonplanar  lifting  surfaces.  The  blade  thickness  plays 
a dual  role,  influencing  both  the  lifting  and  non-lifting  characteristics 
of  the  blade.  As  part  of  the  nonplanar  surface,  it  will  induce  a con- 
tinuous component  of  velocity  on  points  of  the  surface  itself,  thus 
affecting  the  blade  loading  distribution.  On  the  other  hand,  the  symmetri- 
cal flow  disturbance  caused  by  the  blade  thickness  will  influence  the 
pressure  distributions  on  the  suction  and  pressure  sides  of  the  blade 
without  contributing  to  the  net  blade  loading.  Both  effects,  designated 
for  brevity  as  due  to  "nonplanar"  thickness  and  "flow  distortion"  thick- 
ness, will  occur  in  steady-state  flow  conditions  since  the  blade  thickness 
is  a time- independent  quantity. 

The  study  of  the  blade  pressure  distribution  arising  from  both  load- 
ing and  thickness  effects  has  been  undertaken  so  that  the  necessary 
information  can  be  obtained  for  cavitation  and  blade  stress  analyses. 

The  study  treats  both  design  and  off-design  conditions. 
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The  surface  integral  equation  which  relates  the  unknown  loading  to 
the  known  velocity  distribution  on  the  blade  is  solved  by  the  mode  approach 
in  conjunction  with  the  "lift  operator"  technique.  The  selection  of  the 
proper  chordwise  modes  in  the  steady-state  flow  condition  at  design  advance 
ratio  J is  dictated  by  the  shape  of  the  loading  distribution  in  two- 
dimensional  flow  on  a foil  with  the  same  camber  distribution  (NACA  - a mean 
lines,  NACA  mean  lines  of  the  k-  and  5-digit  wing  series,  lenticular  mean 
line,  etc.).  At  off-design  J in  the  steady-state  condition,  there  is  an 
angle  of  attack  due  to  the  difference  AJ  and  the  additional  loading  (to 
that  at  design  J)  due  to  this  angle  of  attack  is  represented  by  the  first 
(cotangent)  term  of  the  known  Birnbaum  chordwise  modes.  The  complete 
Birnbaum  modes,  which  have  the  proper  leading  edge  singularity  and  satisfy 
the  Kutta  condition  at  the  trailing  edge,  are  used  to  represent  the  chord- 
wise  loading  distribution  in  the  unsteady  flow  conditions  at  both  design 
and  off-design  J. 

The  linearized  unsteady  lifting  surface  theory  requires  the  singular 
behavior  at  the  leading  edge.  Although  the  singularity  is  integrable,  its 
presence  in  the  blade  pressure  distributions  is  unrealistic.  Therefore, 
whenever  the  Birnbaum  modes  are  used,  a correction  factor  to  remove  the 
leading  edge  singularity  based  on  Van  Dyke's4  and  Lighthill's5  method  is 
introduced  in  the  blade  pressure  distribution. 

This  research  is  sponsored  by  the  Navai  Sea  Systems  Command  General 
Hydrodynamics  Research  Program,  under  Contract  N0001 A-75-C-0^82 , adminis- 
tered by  the  David  Taylor  Naval  Ship  Research  and  Development  Center. 
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I.  All  OUTLINE  OF  THE  ANALYSIS 

In  the  analysis  which  follows,  the  blades  of  any  ship  propeller  are 
treated  as  warped  lifting  surfaces  which  encounter  the  spatially  varying 
inflow  and,  hence,  develop  unsteady  forces  which  are  cyclic  functions  of 
blade  position.  The  blades  need  not  be  of  a pure  helicoidal  surface,  i.e., 
they  may  have  varying  pitch  over  the  radius.  However,  the  reference  sur- 
face, along  which  the  shed  vorticity  is  considered  to  be  convected,  is  re- 
stricted to  be  a pure  helicoid  whose  local  pitch  is  fixed  by  the  joint 
action  of  the  forward  speed  of  the  ship  and  the  tangential  velocity  at 
any  radius,  or  simply  the  pitch  angle  of  this  vortex-wake  surface  is  tan  ' • 

The  theory  is  formulated  to  give  the  linearized  pressure  distribution 
on  each  side  of  the  blade  at  eight  radial  locations.  These  distributions 
are  made  up  of  an  antisymmetric  part  associated  with  the  lifting  action 
of  the  blade  and  a symmetric  part  developed  by  the  blade  thickness  distribu- 
tion. As  the  blade  surface  is  non-planar,  the  flow  produced  at  any  element 
by  all  other  thickness  elements  produces  a weak  normal  flow  and  this  flow 
is  balanced  by  the  induction  of  ant i -symmetr i c loading  elements  which,  in 
turn,  contribute  to  the  pressure  distribution  as  well  as  to  the  forces 
and  moments.  This  small  loading  is  referred  to  as  being  the  result  of 
the  presence  of  non-planar  thickness;  the  major  part  of  the  asymmetrical 
loading  is  produced  by  the  presence  of  camt  , geometric  flow  angle  and 
the  normal  velocity  components  associated  with  the  hull  spatially  variable 
wake.  The  anti-symmetric  part  of  the  blade  pressure  distributions  are 
solely  responsible  for  the  forces,  whereas  the  symmetrical  and  asymmetrical 
contributions  determine  the  pressures  on  each  side,  the  knowledge  of  which 
is  essential  for  prediction  of  cavitation  inception. 

A.  Propeller  Loading  Distribution 

1.  Integral  Equation 

The  linearized  unsteady  lifting  surface  theory  for  a marine  propel- 
ler, with  its  blades  lying  on  a helicoidal  surface  and  operating  in  non- 
uniform  flow  of  an  incompressible,  ideal  fluid,  is  formulated  by  means  of 
the  acceleration  potential  method.  It  is  based  on  a small  perturbation 
approximat ion  and,  also,  on  the  assumptions  that  the  propeller  blades  are 
thin  and  operate  without  cavitation  and  flow  separation. 
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It  is  known  that  the  pressure  field  generated  by  a lifting  sur- 
face S is  given  by  distributed  doublets  with  axis  parallel  to  the  local 
normal  and  with  strength  equal  to  the  pressure  jump  across  the  surface. 

If  the  acceleration  potential  function,  a scalar  function  of  position, 
is  introduced  and  defined  as  having  its  gradient  equal  to  the  accelera- 
tion vector,  then  the  strength  of  the  doublet  distribution  is  propor- 
tional to  the  discontinuity  of  y between  the  values  for  the  positive 
and  negative  oriented  surfaces  with  respect  to  the  direction  of  the  nor- 
mal n , i . e . , 

W = - V = 7-(P  - PJ  =7-  AP  (1) 

Here  use  is  made  of  the  linearized  relation  between  the  acceleration  po- 
tential and  the  perturbation  pressure,  viz., 


where  P is  local  pressure  and  pf  is  fluid  mass  density.  Thus,  the  pres- 
sure jump  is  defined  as  that  between  the  pressure  at  the  back  (suction 
side)  and  the  pressure  at  the  face  (pressure  side)  of  the  blade. 


The  pressure  P at  any  point  (x,r,b)  (in  cylindrical  coordinates) 
at  time  t due  to  all  pressure  doublets  distributed  over  dummy  (loading) 
points  (•', •>»’->)  will  be  given  by 


P(x,r 


t)  = T~  //  AP  (6 ,0 ,9 ; t) 
S 


J dS 

3n  R'  (x,r,<j>;£,p,0) 


(3) 


where  j_  is  the  normal  derivative  on  the  surface  S at  the  loading 

J point  (",  ,c,G),  with  n the  unit  normal  vector  having  positive 
axial  component 

1 /2 

and  R' (x, r , ; ; ; , p , 0)  = [ (x-6)2  + r2  + p2  - 2rp  cos(9-<j>)]  is  the  Descartes 
distance  between  the  given  control  point  (x,r,<f)  and  the  loading 
point  ( j , p , 0 ) . 

The  relation  between  the  velocity  potential  function  $ and  the 
acceleration  potential  function  has  been  established  by  solving  the 
[ Euler  equation  of  motion: 
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$(x,r,$;t)  = y - / *(t,  r,  <f ; t 


X-T 

U 


■)dr 


00 


where  U is  the  undisturbed  forward  velocity  of  the  lifting  surface  and 
the  integral  represents  the  total  effect  at  time  t of  all  previous  ac- 
celerations of  a fluid  particle  at  present  at  the  observation  point  x. 

For  doublets  with  pulsating  strength  L?  U ,p  , 0)  e ' wt  at  loading 
points  (C,p,0)  which,  together  with  the  observation  (control)  points 
(x,r,q)  are  rotating  with  angular  velocity  -f!  (where  the  negative  sign 
is  introduced  to  accord  with  a right-handed  coordinate  system),  Eq.  (3) 
becomes 


»(x,r,*o;t)  " + ^ //  tPU.fi, Qje  0n  R.  (x , ,p  , 

(5) 

where  w = frequency  in  radians/sec. 

= angular  position  of  loading,  control  points  with  respect 
to  blade  reference  line  in  moving  coordinate  system  (Fig.  2) 

Substituting  Eq.  (5)  in  Eq.  (b)  and  identifying  the  lifting  sur- 

fact  as  the  helicoidal  surface,  x = 6 /a  or  i = 6 /a,  of  an  N-blaced  pro- 

o o 

pel ler  yields 


4>(x,r,$  ;t)  = + 


N 


o'  bnpfU 


iqflt 
r n=l  S 


I e 


ff  AP(q) (E.P.O  ) 


/ e^I^t-x)  -0n]  3 (l)dTdS 
J 3n  R 


(6) 


where  a = ft/U 

q = u/n  = order  of  blade  harmonic 
§ = 2ir.(n- 1 ) /N , n = 1,2,...N 


. 1/2 


R = { (t-e)2  + r2  + p2  - 2rp  cos[eQ  - <t>Q  + Qn  - a (x-x)  ] } 

Then  the  induced  velocity  at  the  control  point  (x,r,$  ;t)  due  to 

o 

the  loading  at  U,p,8  ;t),  when  both  points  are  located  on  the  propeller 
blades,  can  be  determined  from  Eq.  (6).  The  normal  component  of  the 
velocity  induced  by  all  elements  will  be  given  by 


3^r  *U,r, *0;t) 


(7) 
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where  n‘  is  the  unit  normal  vector  on  the  helicoidal  surface  at  the  control 
point,  having  a positive  axial  component. 

The  requirement  for  an  impermeable  boundary  (i.e.,  no  flow  through 
the  lifting  surface)  is 

W + TTT  4>  (x,  r ,4>  ; t)  = 0 

o n o 

where  W is  the  known  ( imposed) perturbat ion  velocity  normal  to  any  blade 
chord 

W - - “ - 4>  (x, r ,<*>  ; t) 

dn  o 


which  expresses  the  equality  of  the  negative  of  se,l  f- induced  normal  velo- 
city with  the  imposed  velocity  in  the  downwash  (positive  n')  direction. 
This  yields  the  integral  aquation  for  the  unknown  pressure  change  AP  at 
frequency  q: 


-e 


iqOt  N 


“--W7T  I //iP(q)(C.P,9„)  ~T  / 
» n-  I j 

where 


o'  3n' 


i q [a (x-x) -0  ] 3 ,1 


£ ^dS 


(3) 


2_  = E (a  J L JL) 

3n  n , 2”  3£  2 

/I  + o 


3 __  r / 9 1 3 \ 

an'  ” yrr r — jrr  d*  ” 2 3$ 

/]  + a^r^  r^  o 


(9) 


It  should  be  noted  that  the  left  side  of  Equation  (8)  is  the  com- 
ponent of  the  flow  normal  to  the  nose  tail  line  of  any  section  of  the 
propeller  blade,  whereas  the  right  side  is  the  induced  velocity  normal 
to  the  helicoidal  reference  surface.  This  lack  of  precise  identity  in 
resolution  of  the  onset  flow  and  the  flow  induced  by  all  loading  elements 
is  in  keeping  with  the  linearized  smal 1 -perturbat ion  theory  employed. 
Specifically,  the  onset  flow  normal  to  the  nose  tail  line  is  (upon  omit- 
ting camber  for  the  present)  expressed  by 


6 
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, u + V 

w - - A T+  )2"+  (fir+VT')2sin(ep  - tan  v~) ) 

where  and  are  longitudinal  and  tangential  perturbation  velocities. 
When  this  is  expanded  to  secure  the  first  order  approximation  (defined  for 
the  condition  ■ 1/7  **  0),  we  obtain 


W ■ - /u1  + (fir)2  s in  (8  -3)  + (V, cosQ  - V^sinG  ) 

P L p T p 


and,  for  small  values  of  S^-S,  this  reduces  to  the  linear  approximation 


W - - »/U2~+  (fir ) 2 (6  -S)  + (V,  cos0  -VT sinO  ) 
P u P 1 P 


Here  S ■ tan  1 = tan  ' — and  6 is  the  pitch  angle  of  the  propeller, 

fir  ar  p 

It  is  to  be  noted  also  that  blades  with  zero  camber,  zero  thickness 
and  having  pitch  angle  variation  6^  = 3 when  set  into  a flow  with  zero 
perturbation,  i.e.,  V^  *»  V_  = O.will  experience  no  crossflow  and,  hence, 
will  develop  no  pressure  loading.  Thus,  all  these  quantities,  incidence 
angle,  camber,  hull  wake  and  inductions  from  non-planar  thickness  are 
considered  as  perturbations  of  the  reference  flow  whose  direction  is  de- 
fined by  8 and  which,  by  itself,  produces  zero  loading.  The  effects  of 
each  of  these  imposed  flows  on  the  blades  are  calculated  separately  and 
simply  added  together  as  allowed  by  the  linearity  of  the  theory. 

The  left-hand  side  of  Eq.  (8)  can  be  written  in  the  form 


W - l V(q) (r,*Q)e  iqQ 


q=o 


(10) 


where  0 = - fit  in  terms  of  a moving  coordinate  system  and  V is  in 
complex  form.  All  but  the  hull  wake  are  perturbat ions  from  the  steacy- 
state  flow  only,  q = 0,  and  are  analytically  computed  from  the  propeller 
geometry.  The  normal  wake  velocities  are  derived  from  an  harmonic  analysis 
of  the  wake  measurements  as  will  be  shown  in  a later  section. 

*See  J.P.  Breslin,  "Determination  of  the  Normal  Velocity  to  a Propeller 

Blade  Section  in  a Ship  Wake,"  Tech.  Memo.  S IT-DL-76-172,  January  1976. 
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On  utilizing  Eq.  (10),  the  integral  equation  (8)  requires  at  each 
frequency  q that 

ap(‘))(£,o,6cj)  K(r,*o;p,Oo;q)dS  (11) 


where  the  kernel  function  is  given  by 

K(r.V»..oi,)-Bi?rl..  I / e'qd (T-x( 

f o->o  n=  1 -«> 

n 02) 

0 ~ 9 

and  6 = — - (x-£)  • 

d 

The  time  factor  has  been  eliminated  from  both  sides  of  Eq.  (11)  and  it 
is  understood  that  the  real  part  of  the  solution  must  be  taken  at  the 
end.  The  limiting  process  is  introduced  in  the  kernel  function  to  avoid 
the  mathematical  difficulty  due  to  the  presence  of  a high-order  singular- 
ity. The  kerne)  function  is  one  of  the  most  complex  of  lifting  surface 
theory  since,  in  addition  to  its  high-order  singularity,  it  presents 
other  complications  arising  from  the  helicoidal  surface  and  from  the 
interference  of  the  other  blades.  Thus,  attention  is  given  to  the  numer- 
ical solution  of  the  integral  equation  by  means  of  a high-speed  digital 
computer.  The  analysis  has  been  carried  to  the  stage  where  laborious 
computations  can  be  efficiently  performed  by  the  numerical  procedure. 

It  has  been  assumed  that  the  shape  of  the  chordwise  loading  dis- 
tribution is  the  same  as  that  in  two-dimensional  flow;  the  spanwise  dis- 
tribution is  left  to  be  determined  by  the  solution  of  the  integral  equa- 
tion. A method  called  the  "generalized  lift  operator"  techn i que1  is  ap- 
plied which  reduces  the  surface  integral  equation  to  a line  integral  equa- 
tion along  the  propeller  radius.  Then,  by  the  collocation  method,  the 
line  integral  equation  is  reduced  to  a set  of  algebraic  equations  in  the 
unknown  spanwise  loading  distribution. 
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The  mathematical  complexity  of  the  problem  has  imposed  a series 
of  concessions  with  regard  to  propeller  geometry  and  helicoidal  blade 
wake.  To  save  time  and  expense  in  computing,  a most  important  concession 
was  that  of  approximating  the  blade  helicoidal  wake  in  "staircase"  fashion,2 
which  appears  to  be  physically  realistic.  However,  with  improvement  of 
the  high-speed  digital  computers,  the  exact  treatment  of  the  blade  wake3 
is  being  uti 1 ized. 

Following  Ref.  (3),  the  substitutions  listed  below  are  made  in 

Eq.  (II), 

L(q)(p,oa)  =Ap(q)  U,P,eo)  • Pe£  (ib/ft)  (13) 


0 = ap  - 0^  cos0 

o b a 

r r 

<i>  = a - 0,  cos$ 

o b a 


0 5 6 < it 
a 


0 1 <J>  $ it 

a 


1-1  Y 

R TT  L 


im[0  -d  +0  -a(t-x)3  00  . / fU 

' °on  / lm(|k|»)Km(ik|r)e'(T-5)k, 


dk  (16) 


for  p < r.  (If  p > r,  p and  r are  interchanged  in  the  modified  Bessel 

functions  I ( ) and  K ( ) of  the  first  and  second  kinds.)  Here  a is 
m m 

the  angular  position  of  the  midchord  line  from  the  generator  line  through 
the  hub  in  the  projected  propeller  plane,  0^  is  the  subtended  angle  of  the 
projected  blade  semichord,  0^  and  $ are  angular  chordwise  locations  of  the 
loading  and  control  points,  respectively,  and  the  superscripts  p and  r 
refer  to  values  at  loading  and  control  points,  respectively. 


After  the  appropriate  chordwise  mode  shapes  Q (n)  are  introduced, 


so  that 


*“*  1 
n = 1 
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where  L^q,n^(p)  are  the  spanwise  loading  components  to  be  determined  by 
the  solution  of  the  integral  equations,  the  chordwise  integration  of 
the  kernel  (Eq.  (12))  over  0^  is  performed.  Then  both  sides  of  the  in- 
tegral equation  (11)  are  operated  upon  by  the  “generalized"  lift  opera- 
tor ; thus, 

l n 

— I $(m)  {Eq.  (11)  d<J>^  (18) 

o 

where 

$(1)  = 1 - (the  Glauert  lift  operator) 

♦ (2)  = 1+2  cos$ 

a 

<t>  (m)  = cos(m-l)4>a  for  m > 2 

A set  of  line  integral  equations  is  thus  obtained  with  maximum 
order  m equal  to  maximum  n of  chordwise  modes: 


■ " r / ,>">(„)  K^hr.p-.qldp 

U n * 1 p 


u<-’”>(r-)-  - i / Hi)  e-i<K>r,  <”)  (q6r) 


/ — \ r , ir  i qO  cos* 

I m (qQ  b)  - — f <t(m)  e ad4>c 


l'  ‘ (y)  is  given  for  various  m in  Appendix  A. 

The  modified  kernel  function  after  the  9^  and  <f>a  chordwise  inte- 


grations is  finally  (see  Ref.  3): 


K(™’'n)  = { ^ } —=—  l e 

i(7rpfU2r  a/l+a2rz  m=-°° 

• o i 


- imAa 


m=q+£N 


• {eiMAa[a2(m-q)  +-2-]  [a2(m.q)+JL]  | ^ (qef)  A( ("}  (qO^) 
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I 

I 

I 

I 

I 

I 

I 

I 

I 

I 


' I (alm-qlp)  K (alm-qlr) 
m m 


7 / !“  ^^5 — '.(IMklK^lklDdk) 

-«  rz  pz  v ^ 


(n) 


where  A (z)  = / J (n)e 
o 


-iz  cos0 


(21) 


asir>0  do  and  Ac  = cjr-a°. 
a a 


In  Eq.  (21),  all  terms  outside  the  first  braces  are  nond imens ion- 

alized  with  respect  to  propeller  radius  r as  is  also  p in  Eq.  (19). 

/ ” \ _ o 

A ' 1 (z)  for  var ious  n and  0 I'n)  is  also  given  i n Appendix  A . 

The  integral  equations  given  by  Eq.  (19)  with  the  kernel  given 
by  Eq.  (21)  are  solved  numerically  by  the  usual  collocation  method,  with 
the  loading  L^q,n^(p)  assumed  to  be  constant  over  each  small  radial  strip. 
Then  only  the  kernel  needs  to  be  integrated  over  the  radial  strip. 


Reference  3 details  the  analytical  development  and  the  various 
numerical  procedures  to  obtain  the  finite  contributions  of  the  Cauchy- 
type  s i ngu  1 ar  i ty  of  the  k-integral  at  k = a(m-q)  and  of  the  higher  order 
Hadamard-type  singularity  when  p * r. 

2.  Chordwise  t-lodes 

The  proper  selection  of  chordwise  modes  is  dictated  by  the  loading 
distribution  on  a foil  in  two-dimensional  flow.  Thus,  in  the  unsteady 
flow  case,  the  unknown  loading  function  is  approximated  in  the  chordwise 
direction  by  the  known  Birnbaum  distribution  which  has  the  proper  lead- 
ing edge  singularity  and  satisfies  the  Kutta  condition  at  the  trailing 
edge.  In  this  case,  Eq.  (17)  is 


/ \ n max  / ~ \ 

L(q) (P,6a)  = J L(q’n)(p)  0 (n) 

n = 1 

, / i \ 9 n max  , -» 

= a(q,1)  cot-f  + l L(q’n)(p)  sin(n-l)9  } (22a) 

1 n = 2 a 


In  the  steady-state  flow  condition,  the  chordwise  mode  shapes 
are  selected  to  conform  to  the  observed  pressure  distributions  of  the 
NACA  foil  sections.  Furthermore,  the  analysis  and  corresponding  program 
is  divided  into  two  parts,  for  design  and  off-design  propeller  conditions. 


I 

I 


11 


1 


R-1869 


At  design  advance  ratio  J,  the  propeller  operates  at  almost  optimum  ef- 
ficiency and  the  cbordwise  loading  distribution  due  to  blade  camber, 
thickness,  angle  of  attack  and  inflow  conditions  (mean  wake  effects) 
does  not  deviate  much  in  shape  from  the  loading  distribution  of  an  air- 
foil in  two-dimensional  flow.  At  off-design  advance  ratio,  however, 
the  same  propeller  in  the  same  wake  is  subjected  to  a change  in  angle 
of  attack  due  to  AJ , the  difference  between  off-design  and  design  ad- 
vance ratio.  The  additional  load  distribution,  due  to  this  change  in 
angle  of  attack  in  the  off-design  condition,  is  in  accordance  with  thin 
wing  theory  represented  by  the  first  term  of  Eq.  (22a). 


The  camber  distribution  is  a decisive  factor  in  the  selection 
of  the  proper  mode  shapes  of  the  chordwise  loading  distribution  in  the 
steady-state  flow  condition,  since  it  is  known  that  many  properties  of 
wing  sections  are  primarily  functions  of  the  shape  of  the  mean  line. 

Thus,  the  program  in  the  steady-state  case  deals  with  two  different  types 
of  mode  shapes. 

The  first  is  appropriate  to  the  NACA-a  mean  lines  characterized 
by  the  fact  that  the  chordwise  loading  is  constant  from  x/c  = 0 at  the 
leading  edge  to  x/c  = a (a  varying  from  0 to  1 ) and  then  decreases  linear- 
ly to  zero  at  x/c  = 1,  the  trailing  edge.  For  brevity,  this  type  is 
designated  as  "roof-top"  loading.  In  this  case,  Eq.  (17)  becomes 


L(o)(p.0a)=L(p)(°’,)  0 (1) 


-e 


L(p)  (o,1).  0 S x 5 a 

(p)(o,1)  (~)  , a f x i 1 
l -a 


(22b) 


and  A ^ (z)  of  Eq.  (21)  is  as  gi ven  in  Append i x A 


The  second  type  of  mode  shape  in  the  steady  flow  case  is  appro- 
priate to  NACA  mean  lines  of  the  !>-  and  5“digit  wing  series  (including 
the  lenticular  mean  line)  and  to  arbitrary  mean  lines  in  general  which 
are  characterized  by  the  fact  that  the  chordwise  loading  is  finite  every- 
where and  zero  at  both  edges  and  can  thus  be  described  by  a series  of 
sine  mode  shapes: 


. (o)  , . \ I n nlaX  . / ■,  (o,n)  . - 

L (p’ V “ 7 l L(p)  sin  n ea 


n = 1 


(22c) 
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It  is  to  be  noted  that,  when  the  complete  Birnbaum  mode  series,  Eq.  (22a) , 
or  the  cotangent  term  alone,  is  used  and  the  solutions  for  the  span- 
wise  loadings  are  substituted  in  the  chordwise  distribution,  a correction 
for  the  leading-edge  singularity  is  introduced  based  on  Van  Dyke's  and 
Lighthill's  method4*5. 

3.  Perturbation  Velocity  Distributions 

As  mentioned  earlier,  the  left-hand  side  of  the  integral  equa- 
tion represents  the  positive  normal  components  of  the  velocity  perturba- 
tions from  the  basic  flow  due  to  non-uniformity  of  the  flow  (wake),  to 
blade  camber,  blade  thickness  (non-planar)  and  incident  flow  angle,  in 
the  design  condition,  and  to  an  additional  angle  of  attack  in  the  off- 
design  condition  arising  from  the  difference  between  off-design  and 
design  advance  ratios. 

a)  Normal  Velocity  Due  to  Wake 

This  is  the  contribution  of  the  flow  non-uniformity  to  V^(r)/U 
of  Eq.  (20).  The  positive  normal  component  (i.e.  downwash)  of  the  wake 
velocity  along  the  middle  chord  of  each  radial  strip  is  given  by 

VM(r)  VL(r)  V (r) 

___  = __  cosep(r)  . sin9p(r) 

where 


V,  (r)  V 

-L = -'  + TT 

U U 


perturbation  of  ur.,form  flow 


Vjr) 


= measured  axial  velocity 


VT(r) 

0 


= measured  tangential  velocity 


0 (r)  = tan"' 
p 2*rro 

P(r)  = propeller  pitch,  ft. 
rQ  = propeller  radius,  ft. 

r = radial  distance  in  terms  of  r 
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so  that 

V)  VJr)  VT(r) 

— jj [— u~  - •]  cos0p (r)  ij—  si nep(r)  (23) 

It  is  to  be  noted  that  the  present  sign  convention  requires  that, 
looking  forward,  the  positive  position  angles  0 are  defined  counter- 
clockwise from  0 * 0 at  the  upright  position  of  the  blade,  that  the 
axial  component  of  the  wake  is  positive  downstream  (positive  x-direction) 
and  the  tangential  component  Vy  is  positive  in  the  counterclockwise  direc- 
tion. 

The  harmonic  analysis  of  the  measured  wakes  yields 


Vr) 


a ( r)  + l [a  (r)cosq  0 + b (r)sinqe] 
° q=l  Q q 


Vy (r)  « 

— jj — = AQ(r)  + l [A  (r)cosqe  + B (r)sinq  Q ] 
q=l  q s 

The  zero  harmonic  of  the  normal  velocity  is  then 


v'0)(r) 


(a  (r)  - 1 ) cos0  (r)-A  (r)  sinO  (r) 


and  the  other  harmonics  are 


V (r) 

— — -jj = [aq(r)  cos9p(r)  - Aq(r)  s in9p (r)]  cosq9 


+ [bq(r)  cos9p(r)  - Bq(r)  s in9p (r)] s i nq9  (23b) 


As  stated  previously  (see  Eq.  (10),  the  left-hand  side  of  the 
integral  equation  must  be  in  complex  form 


[(~)o+i(J^.]e'iq0 


U 'R  U 
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and  it  is  understood  that  the  real  part  is  to  be  taken  eventually.  The 
real  part  is 

w (q)  v(q) 

* N N 

( — u — ) R cos  q 0 + (— J- ) , S i n qo 

Therefore,  the  input  to  the  program  for  the  real  part  is  from 
Eq.  (23b) 

(r) 

( rr ) = a (r)  cos0  (r)  - A (r)  sinO  (r)  (24) 

u K q p q p 

and  for  the  imaginary  part 

v'^fr) 

( rr ),  = b (r)  cos9  (r)  - B (r)  sine  (r) 

J I q p q p 


b)  Normal  Velocity  Due  to  Blade  Camber 

The  velocity  normal  to  the  blade  in  negative  direction  induced 
by  the  flow  disturbance  caused  by  blade  camber  is  independent  of  time 
(since  the  blades  are  considered  rigid)  so  that  only  the  steady  state 
loading  is  affected.  The  loading  due  to  camber  effect  is  obtained, 
therefore,  from  the  steady-state  part  of  integral  equation  (19)  with 
V^°  (r)/U  of  (20)  replaced  by 


(r) 


U 


/nw  ayr,s) 

c(r)  3s 


_ a— — 5 — r 3f  (r.d,  ) 

2m  + _£__o_ 

c(r)  sin  <{i  34> 

a a 


where 


c (*"  ,s) 

j slope  of  the  camberline  fc(r,s)  given  at  discrete 

55  points  measured  from  the  face  pitch  line 

s = (I  - cos<!>  )/2,  chordwise  location  non-d  imens  iona  1 ized 
on  the  basis  of  c(r) 

c(r)  = chord  in  feet 

4>  = angular  chordwise  position 


|' 
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On  application  of  the  generalized  lift  operator  (see  Eq.  ( 1 8 ) ) 
to  both  sides  of  the  integral  equation,  the  left-hand  side  of  Eq.  (19) 
for  this  particular  component,  camber,  becomes  for  each  order  m 


y (o,m) 
c 


(r) 


U 


2/l  + a^r- 

~7UV~ 


TT 

/ $(m) 
o 


Of  d<p 

C 'i 

06  sin^ 
a a 


(26) 


where,  as  can  be  seen  from  Eq.  (18),  4>(m)  can  be  expressed  in  terms 
of  cosm^.  Integrating  by  parts 

1 , <*f  f (r,s)cosmb 

I r cosnvj) c , _ c 

c(r)  ‘ sing  3g  c7r)sing 


+ 


fc(r,s) 
s in  2g 


[m  sinmg  s i n<J>  + cosm<jj  cosg]dg 


(27) 


In  the  small  range  near  the  leading  edges,  0 ~ - cos  '(0.9)  or  0-  si  0.05, 

the  camberl  ine  is  assumed  to  be  parabolic  and,  in  the  range  cos  ' (~  .S)^  <j>  5 tt 
or  0.90  t s i 1.0  near  the  trailing  edge,  the  camberline  is  assumed  to 
be  a straight  line.  The  integration  is  done  analytically  in  these  re- 
gions and  numerically  over  the  remainder  of  the  chord.6  The  required 
input  information  is  the  ratio  of  camber  ordinate  fc(r,s)  to  chord  c(r) 
at  19  equidistant  positions  from  s = 0.05  to  0.95  of  chord.  For  arbi- 
trary camber,  these  ratios  are  read  into  the  program.  The  program  has 
options  to  compute  them  in  the  cases  of  NACa-a  meanlines  (a =0.8,  0.8 
modified,  1.0)  and  of  lenticular  (s  ine-sauared)  camberline  from  the  ratio 
of  maximum  camber  ordinate  to  chord  at  each  radial  position,  and  in  the 
case  of  the  NACA-l*  digit  series  from  the  ratio  of  maximum  camber  ordinates 
to  chord  and  the  chordwise  position  of  the  maximum  ordinate. 


o)  Normal  Velocity  Due  to  "non-planar"  Blade  Thickness 


The  thickness  of  propeller  blades  describing  a helicoidal  or  non- 
planar  surface,  as  stated  earlier,  generates  a velocity  field  on  the 
blade  itself  and,  consequently,  affects  its  loading  distribution  (in 
the  steady-state  only  since  the  blades  are  considered  rigid).  The  nor- 
mal velocity  component  due  to  the  effect  is 


VT(o)  _ , 3$T(x,r,*o;0) 

”u~  " u aTT* 


(28) 
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where  3/3n‘  is  as  defined  in  Eq.  (9). 

Resorting  to  the  "thin  body"  approximat ion,  the  velocity  poten- 
tial <J>t  is  given  by  that  of  the  source-sink  distribution  with  strength 
proportional  to  the  slope  of  the  thickness  distribution.  In  the  case 
of  an  N-blade  propeller,  the  velocity  potential  is  given  by 
1 N M(€,p,0  ) ,fy~T — 2 — T 

--ill  I ■ - °-  ‘>dpd()o  <« 

n=!  0 p 
o 

Here  the  source  strength  M is 


M(?,p,0o)  = 2U 


df  (p,o)  3f  (C,P,e  ) 

^-TT— 2U  -Hr-- 


where  f (?,p ,0  ) = thickness  distribution  over  one  side  of  the  blade 

T O 

section  at  radial  position  p in  the  propeller  plane. 
The  Descartes  distance  R is 

R = { (x-£)2  + r2  + p2  - 2rp  cos  (0q  - + 6^) } 

where  x = /a  = (or  - 0["cosc}>  )/a  , 0 f ^ S 

o b a a 

E = D /a  = (a'L  oPcose  )/a  , 0 < 0 < - 

o b a oi 

0n  = 27r(n-l)/N  , n - 1 ,2,. . .N 

If  the  reciprocal  of  R is  again  expressed  in  a series  as  in  Eq.  ( 16), 


-i— (!)  = i 

3n'V  it 


00  im(0  -<p  +6  ) ■ f r\< 

1 e 0 ° " / (ak+-2-)(lK)m  e'(x'C)kd 

is:  — ao  — co  K 2 


where  (IK.)  = f I ( | k | p ) K (|kjr)  , p < r 

m J m m 

1 >m(lkl r)Km(|k|p)  , P > r 

N im0  fN  for  m = J.N,  1 = 0,  -1,-2,... 
Since  l e n =< 

n=!  I 0 for  al 1 other  m 


and  de  =oPs  i nO  d6  and  —r=  


o b a 


3E  QP  • a 39 
0,sm0  a 
b a 


Equation  (28)  can  be  reduced  to 
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V<o)  (r) 


r -imAo  "r  t Df  (p,0  ) 
l e 11  — L_ — /T+  "a^n2 


. / (ak  + ~)  ( , K)  eiA°l</aei(k-am)0Pbcos0a/a 


. e"«  (k-am)0[cos<Pa/a 


dk  dij  d9 


Before  the  0 - i n tegra t ion  can  be  performed,  it  is  necessary  to 

supply  an  analytical  expression  for  the  thickness  distribution  f (p,0  ). 

t a 

When  the  NACA  k and  5 i g 1 1 wing  series  were  derived,  it  was  found  that 

the  ratio  of  thickness  distribution  y to  maximum  thickness  t of  effi- 

t o 

cient  wing  sections  were  nearly  the  same  when  their  camber  was  removed 
(mean  line  straightened).7  The  thickness  distribution  of  these  sections 
is  given  by 

y t 

— 0 - • 630x  - 1 .758x2  + 1 .A215x3  - .5075x4) 

where  c is  the  chord  length,  x is  the  chordwise  location,  0 1 x f 1,  and 

the  coefficient  of  /x  is  equal  to  /2p  , p being  the  ratio  of  the  lead- 

o o 

ing  edge  radius  to  the  chord  c.  In  general,  for  any  propeller  the  blade 
thickness  distribution  in  the  projected  plane  can  be  approximated  by 


t (p) 

fT  (p,ea)  ~ 2p0P  {/2pQ(p)  /x+  c(~]~  [ax  + bx2  + cx3  + dx4]  } 


1 “CO 50 

and  with  x = r — — 0 < 6 in 

2 a 


fT(p»ea)  “ 2p6P 


a o 
s,n  T+T1 


t (p) 


-[ao(p)  + [ a (p)cos  n 0 ] } 


where  f^  and  p are  fractions  of  propeller  radius  rQ,  Pq(p)  is  the  ratio 

of  leading  edge  radius  to  chord  length  c(p)  at  a given  radial  position 

and  t /c  is  the  ratio  of  maximum  thickness  to  chord  at  that  radial  posi- 
o 

tion.  The  coefficients  of  the  fourth-order  polynomial  in  x are  calcul- 
ated by  the  “least  squares"  method  and  then  with  the  trigonometric  sub- 
stitution for  x,  the  coefficients  an  of  Eq.  (3')  are  obtained  and  inserted 
in  the  input  to  the  program. 
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The  slopes  of  the  thickness  distribution  are,  from  Eq.  (31) 
df  (p.e  ) 9 „ T (p)  U 

~I^~^=p9Pby2^coSf  -2pebTT?r  l nan(p)sinnea)  (32) 
a n=  l 


Lett i ng 


n 3f  (p , 6 ) i ((k-am)9^/a)cos9 

- / V"e  b 

o a 


de 


(33) 


and  applying  the  generalized  lift  operator,  Equation  (30)  becomes  for 
each  lift  operator  mode  m 


V(o,m)(r) 

T 

U 


iNr 


l 


2it2  /l  +a^r^  m=-°° 
m=£N 


rimA0  / /UaV 


f (ak  + -~)  ( I K)  eiAokyal(m)((am-kjer/a)A((k-am)6f’/a)dkdp 
J ■>  m d o 

r Ma) 


Let  u = k - am,  then 
r,(o,m) 


v;w”"'(r) 
0 


iNr 


27t2/  1 +~azr2  p 


/ /r+'a^ip2' 


£ / [au  + lH(a2  +-)]fjlN(|u  + aLN^)K£N(|u  + atNjr)elAau/a 

= -co  -co  r2 

■ l ^ (-u0r/a)A(u6^/a)dudp  (3^b) 

0 0 


(for  p 


U 


<r). 

Finally 

ilL *1 / /TTFp2’  / lm.Part{eiWa!(^(-uS:/a)A(u9D,/a)} 

p o b b 

• {au I (up)K  (ur)  + £ { (au  + tN(a2  + — ) ] ( N ( |u  + a£N  j p) K N ( | u + a£N|  r) 
O O jjij)  r2  x,r‘ 

+ [au-£N(a2 +■—)]) lN(Su-aHN | p)K^N(Su-a£Njr)))dudp  (3^c) 


i 

J 
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{for  P<r,  otherwise  p and  r are  interchanged  in  the  modified  Bessel 
functions.)  As  can  be  seen  from  Appendix  B 


t (p) 


A (A)  =P0£  h Po(p)  B (X ) -2p0^  {2ax  (p)G(X)  + 3a3  (p)  [6G  (A)  - -y-  F(X)] 


where  B(A)  = 


n 

/ e 
o 


cos  Y 


+ i 8a2(p)  F (A) 


(35) 


32a  (p ) . 

+ i 1 [1.0(A)  + (A  - ~)  F(A)  ] } 


(See  Appendix  B,  5) 


G(A) 
F (A) 


s inA 


sinA  - A cosA 


(It  can  be  shown  that  when  u = 0 the  integrand  of  (34c)  is  zero.) 

d)  Normal  Velocity  Due  to  Flow  Angle 

The  incident  flow  angle  has  been  defined  as  the  difference  between 

the  geometric  pitch  angle  0p(r)  = tan  1 (P(r)/27rr)  and  the  hydrodynamic 

_ ) v 

pitch  angle  B(r)  = tan  (1/ar)  of  the  assumed  helicoidal  surface.  The 
normal  velocity  (in  the  positive  direction)  v|°^(r)/U  due  to  this  effect 
depends  only  on  the  radial  position,  not  on  the  chordwise  <J>a  position, 
and  is  given  by 


Vf0)(r)  

- -AT  a2r2(ep(r)  - 6 (r)  ] 

Therefore,  after  application  of  the  generalized  lift  operator,  it  becomes 


_ 


= f -A  + a^r2  [ 8 (r)  - 6(r 
4 0 


)]  for  m = 1 or  2 
for  m > 2 


(36) 


s i nee  I ^ (0)  = I ^ (o)  = l and  I > ^ (0)  ■ 0.  (See  Append  ix  A ) 


In  Eq.  (36)  a = ir/J  where  J is  design  advance  ratio.  At  off-design 
advance  ratio  the  propeller  is  subjected  to  a change  in  angle  of 

attack  due  to  AJ  = - J.  The  normal  velocity  perturbation  due  to 

this  change  in  angle  of  attack  is 
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« {-/T+aJ.r^e  (r)  - tan"'  + 

od  p it  r 


frT[ep(r)  - B(r)] (U/UQd)} 


for  m * 1 or  2,  where  aQd  = anc*  ^ocj  's  ^e  off'design  forward  velocity. 

4.  Blade  Pressure  Distributions,  Hydrodynamic  Forces  and  Moments 
and  Blade  Bending  Moments 

a)  Pressure  Distribution 

After  the  integral  equation  (19)  is  solved  for  propeller  spanwise 
loading  components  L^q,n^(p)  (lb/ft)  due  to  the  various  perturbation 
velocity  distributions,  the  spanwise  loading  distribution  is  determined 
from  Eq.  (17): 

L(q)  (r)  = / L(q)  (r ,0  )sin8  d6 
' ci  a a 

o 

tt  n max  , -* 

= / A Liq’n)(r)  0(n)sin9ad0a  (38) 

o n * 1 

For  the  unsteady  loading  due  to  wake  (q  t 0) , the  complete  Birnbaum 
series  is  used  with  the  chordwise  loading  distribution  given  by  Eq.  (22a). 
The  spanwise  loading  distribution  is,  in  this  case, 


L(q)  (r)  =-/ (L(q,1)  (r)  (1  +cos0  ) + l L (q  ,n)  (r)  s in  (n-1 ) 6 s i n6  }d0 

TT  _ QL  Ct  u u 


= L(q,,)(r)  l(q,2)(r)  (39a 

For  the  steady  loading  due  to  AJ  in  the  off-design  condition,  the  span- 
wise  loading  distribution  is  the  first  term  of  Eq.  (39a). 

In  the  steady  state  case  (q = 0)  for  NACA  a mean  line  blade  sec- 
tions, the  chordwise  loading  distribution  is  the  so-called  "roof-top" 
distribution,  given  by  Eq.  (22b),  which  can  also  be  expressed  as 
, (o,l)(r)  for  o < 0^  < cos  '(l-2a) 


, (o)  , . < , . \ 1 + cos8  , 

L (r)  =[L(°.1)(r) a fQr  cos_1  (1  -2a)  <e  < * 

2(1  -a') 

The  spanwise  loading  distribution  in  this  case  is 
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■I 


C \ C i \ cos  (l-2a)  , r 

L ° (r)  =L  °*  (r){  / sinO  d9  + ./  (l+cos6  )sinu  dG  } 

i a a „ , . - » - 1 v , ~ - \ ct  a a. 

0 2(l-a)  cos  Cl -2a) 


= L ) 


(r)  (a  + 1 ) 


(39b) 


The  steady  chordwise  loading  distribution  for  NACA  mean  lines 
of  the  4-  and  5"digit  wing  series  and  for  arbitrary  mean  lines  in  gen- 
eral is  given  by  Eq.  (22c).  The  spanwise  loading  distribution  is  then 
given  as 


L(°)(r)=i/J  L{o’R)(r) 


- 2 L 


o n=l 

(o,  1 ) 


s in  n G s i nG  d6 
a a a 


(r) 


(39  c) 


The  chordwise  pressure  distributions  in  lb/ft2  are  obtained  by 
dividing  the  chordwise  loading  distribution  by  the  semi-chord  in  feet 
(see  Eq.  (13)). 


Whenever  the  cotangent  Birnbaum  mode  is  used,  the  Van  Dyke  or 
Lighthill  correct  ion4 > 5 is  applied  to  the  chordwise  pressure  distribu- 
tion to  remove  the  physically  unrealizable  leading  edge  singularity. 


Their  simple  rule  for  obtaining  a uniformly  valid  solution  is 


'Pt  S+Po/2  P> 


where  C0  = pressure  coefficient  according  to  first-order  thin  airfoil  theory 


s = (l-cos0  )/2  = chordwise  location  from  the  leading  edge  in 
a 


terms  of  chord  length 


= leading  edge  radius  of  the  profile  in  terms  of  chord  length 


The  modified  chordwise  loading  distribution  is  then 


L(q)(r,0j  -1  (TTT^~ ^-){L(q,1)(r)cot  f+  l L (q>n)  (r) si n (n-1 ) 0^ 


a t r 1 + p - cos0 

o a 


n=2 


(40) 


Since  the  pressure  near  the  leading  edge  is  governed  by  the  first  term 
of  (40) , viz. , 
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I 


hi 


I 

I 

I 

I 


1 - cosO  0 s inC 

(t  _u_  _ a 

I +0  - cosO  CO  2 1 p - cosO 

o a o ct 


the  location  of  maximum  pressure  near  the  leading  cage  is  determined  appar- 
ently from 

S i n 0 

d q ^ 


dG  1 + p - cosO 
a o a 


Thus  (cosO  ) = 1/0+P  ).  This  is  the  chordwise  location  of 

a max  o 

a point  where  the  maximum  pressure  near  the  leading  edge  will  occur  after 

the  application  of  the  Van  Oyke-L i ghth i 1 1 correction,  provided  pQ  is  not 

too  larqe.  If  p > .015,  the  sine  series  of  (*40)  will  influence  the 
o 

location  of  maximum  pressure  near  the  leading  edge. 


b)  Propeller-generated  Forces  and  Moments 

The  principal  components  of  the  hydrodynamic  forces  and  moments 
which  are  evaluated  are  listed  below  and  shown  in  Figure  1 with  the  sign 
convention  adopted: 


Forces : 


F 

x 

F and  F 

V z 


thrust  (x-d i rect ion) 

horizontal  and  vertical  components, 
respectively,  of  the  bearing  forces 


Momen  t s : Q 

Qy  and  Q.z 


torque  about  the  x-axis 

bending  moments  about  the  y-  and  z-axis, 
respect i vel y 


The  total  forces  at  frequency  UN U = 0, 1 ,2. . . ) induced  by  an  N-bladed 
propeller  are  determined  as  (see  Appendix  C) 


F 

y 


Re  { Nr  eaNfit  / (r)coso  (r)dr) 

° o p 

Re  } I [LUN-,’"5(r)A(n")(-0^)+LUN+,*")(r)A^)(6^] 

2 ° n=l 

sinO  (r)dr) 

P 
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F -Re  {^eiMint  / £ [L(WM  *">  (r)A("}  (-6^)-  L {*N+I ’">  (r) A(">  (e!")] 

z 2 1 ' b b 


n“l 


sin9  ( r ) d r } 
P 


(Ma) 


The  moments  are  determined  by: 

I 

/ 

o 


Q = -Re  {Nr2  e'lNfit  j L^N^(r)sin9  (r)rdr} 
x ° o P 

Nr2  1 - - » 

Q,  =Re{-^eUNflt  / { £ [LUN‘,*n)(r)A(n){-9p+LUN+1’n)(r)A(n)(9^)]cose(r) 

' ° n-1 

+ £ [LUN  ' ,n^  (r)  A^  (-eT)' L'lN+ ' ’ n^(r)  a/R^  (e|")]  ( ie!”)  s inO  (r)tanG  (r)}rdr; 

lb  'bop  p 

n=l 

-Mr-2  i « _ _ 

Q .Re^e1"*11  / { J [L(“'l'")(r)A(,')(-e")-L<'N+l’")(r)»('’)(6")JcoSe  (r) 

2 i i ' '•  b bp 


o - . 
n=  1 


+ £ [L(iN‘,'n)(r)A1(n)(-e^K(J'N+1’n)(r)A1{n)(6[)](ie[)  s in9p  (r)  tan0p  (r)  }rdr, 

n=I  . (4l  b) 

where  A (z)  and  A^  (z)  are  as  defined  in  Appendix  A. 

It  is  seen  from  Eqs.  ( 4 1 a , b)  that  the  transverse  bearing  forces 
and  bending  moments  are  evaluated  from  propeller  loading  components  L^q,n^(r) 
associated  with  wake  harmonics  at  frequencies  adjacent  to  blade  frequency, 
i.e.  at  q = UN  - 1 , whereas  the  thrust  and  torque  are  determined  from  the 
loading  L^(r)  (given  by  Eqs.  39a, b,c)  at  blade  frequency  q = UN.  At  1=  0 
(steady-state)  the  mean  transverse  forces  and  bending  moments  are  deter- 
mined at  shaft  frequency.  Thus: 

Nr  1 


F„-  Re  {-^  /[  £ L(1’n)(r)A(n)(e£)]sin9p(r)dr} 


o - , 
n=l 


Nr  1 


F,«  Re  {-jf  /[  £ L(,'n)(r)A(n)(8[)]sinep(r)dr} 


° *‘-1 
n=  I 


L 
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Nr  1 1 / . ~s  /~\ 

^ / ([  l L (r) A " (0^)  cosGp(r) 

° n=  1 


[ £ L^1  ,n^  (r)A^  (Op]  (iG^)sinOp(r)taaup(r)}rdr; 

n=  1 


NR  2 1 n -x  - 

Q = Re  v-^2-  / {[  £ Lv,,  ,)(r)Al  ' (i  ‘ ) cosO  (r) 
z 2 1 J ? 

° n=l 


[ [ L(l>n)  (r)/i1(n)  (0^)](iGpsir,0p(r)tan0p(r)}rdr;-  (l*l c) 


n=  I 


All  forces  and  moments  can  be  written  in  the  form 


Rg[C^e'  q e 'qQ]  = C ^cos(-qO  + <>q)  = C ^ cos  (qb  - y^) 


where  0 is  blade  angular  position,  positive  in  tiie  cour.ter-clockjise 
direction  from  zero  at  the  upright  position  (12M),  q is  order  of  shaft 
frequency,  is  magnitude  of  force  or  moment,  and  is  phase  angle 

(electrical)  determined  as  the  angle  whose  tangent  is  the  imaginary  part 
(sine  part)  over  the  real  part  (cosine  pari'.  (The  program  output  gives 
and  4>  .) 

q 

In  comparing  phase  angles  with  those  of  experiment,  care  should 
be  taken  to  ascertain  if  measured  values  refer  to  lead  or  lag  position 
and  if  the  quantity  considered  is  at  peak  or  trough.  When  the  theoret- 
ical <i>  is  positive,  the  peak  (or  trough)  of  the  quantity  is  to  the 
left  of  the  upright  position,  i.e.,  leads.  When  is  negative,  the 
peak  (or  trough)  is  to  the  right  of  the  upright  position  and  lags  that 
of  input. 


I 

I 


25 


R-1869 


a)  Blade  Be>iding  Moments 

The  blade  bending  moment  about  the  face  pitch  line  at  any  radius 

r.  of  a blade  is  calculated  from  the  chordwise  integrated  loading  (span- 
J (Q\ 

wise  component)  L M (r)  at  any  shaft  frequency  q as 


iqftt  7 iqftt 
e = r z e 


l^(r)  cos[0  (r)  - 6 (r . ) ] (r-r . )dr 


P J 


(42) 


The  positive  blade  bending  moment  about  the  face-pitch  line  is  that 
which  puts  the  face  of  the  blade  in  compression  (see  Figure  1). 


The  instantaneous  blade  bending  moment  distribution  when  the  pro- 
peller swings  around  its  shaft  in  the  clockwise  direction  is 


M 

‘b 


* I!Mq'jCOs(q0  -4>q) 


(43) 


where  is  the  phase  angle  (electrical)  determined  by 


d 


q 


tan 


(Mq) 

b I m 

(Mq)Re 


and  (M^)ge>(^^) |m  are  rea'  and  imaginary  parts  of 

It  should  be  noted  that  in  the  program  the  value  of  r.  of  Eq.  (k2) 
is  limited  to  any  of  the  midpoints  of  the  radial  strips  into  which  the 
blade  span  is  divided,  at  which  points  the  pitch  angles,  as  well  as  other 
geometrical  characteristics,  are  given  as  input.  The  bending  moment  at 
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any  other  radial  position  can  be  obtained  by  interpolation  or  extrapola- 
tion. For  example,  for  rj  < (r^  - Ar),  where  r;  is  the  midpoint  of  the 

initial  radial  strip  from  the  hub,  r.  , and  Ar  is  the  span  of  the  radial 

n 

strip  (rj  ■ r^  + Ar/2),  the  bending  moment  can  be  estimated  as 


[M  {r  ) -M  (r  = r +Ar)] 

W-  "bV  - ■--■■-V"1 Vr.> 


For  the  case  of  bending  moment  about  the  root  r,,  this  formula  gives 

n 


Hb(rh>  ■ "b<ri>  + 


Hb<ri'-"b<r2» 


5.  Frictional  Thrust  and  Torque 

To  obtain  the  frictional  contribution  to  thrust  and  torque,  one 
must  resort  to  empirical  formulae  for  the  frictional  coefficient  in  tur- 
bulent flow  since,  in  this  case,  the  theoretical  analysis  is  extremely 
difficult.  Hoerner8  suggests  a formula  based  on  statistical  analysis  of 
measured  velocity  distributions  across  the  boundary  layer,  viz.  the  Prandtl- 
Schlichting  formula  for  the  friction  coefficient  for  one  side  of  a smooth 
flat  plate  of  length  c: 


CF  = 


0.455 

(,°9lORc) 


2.58 


(44) 


where  Rc  is  Reynolds  number  based  on  length  c (of  the  expanded  chord) 
i.e.,  in  the  nomenclature  of  this  paper, 

r (r)  = - ZEIjIZI  c(r) 

where  v is  kinematic  viscosity. 

Because  of  thickness,  the  average  velocity  around  a symmetrical 
foil  section  is  higher  than  that  of  the  undisturbed  flow,  the  increase 
being  proportional  to  the  thickness  ratio  tQ/c.  For  the  conventional 
section,  the  total  frictional  drag  coefficient  is  approximated  by8 

t (r) 

CQ(r)  = [2  CF(r)](l  +2  -^-]  (45) 
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Although  no  al lowance  is  made  for  viscous  pressure  drag,  nor  drag 
associated  with  the  deviation  of  the  local  flow  direction  from  that  along 
the  nose-tail  line,  this  formula  agrees  basically  with  theoretical  analyses 
of  skin-friction  drag  and  is  substantiated  by  experimental  results8. 

Taking  AF^(r)  (elemental  friction  force)  as  the  friction  force  per 
foot  of  propeller  blade  span  with 

AFq ( r)  = --  pfU2  (l+a2r2)  c(r)  CD ( r)  (46) 

for  an  N-bladed  propeller,  the  total  force  in  the  positive  x-direction 
and  moment  about  the  x-axis  due  to  skin-friction  drag  are 


F0x  = N / AF0(r)  sine  (r)  dr 
and  r 

QQx  = N / AFD(r)  cosS  (r)  r dr 

r ^ 


On  substitution  of  (44)  - (46)  and  noting  that  the  expanded  chord  length 
is  equivalent  to 


c (r)  = 


2 r r 9, 
o b 

cos's  frT 
P 


where  r is  a fraction  of  propeller  radius  r , the  frictional  thrust  and 
torque  become 

I t (r) 

F0x  = Nro  pfU2  /2Cp(r)t'  +2  c°(r)  '1  0 +a2r2)r0^  tanG  (r)dr 

m 

1 t (r) 

QDx  = Nr^  PfU2  / 2Cp(r)  [1+2  ] (1  +a2r2)r26^  dr 


F.  and  Q.  are  added  algebraically  to  F and  Q , respectively, 

Ox  Ox  ' x x 

as  calculated  in  the  preceding  section.  Thrust  will  be  decreased  slightly 
and  torque  increased  more  noticeably. 

As  was  mentioned  earlier,  this  is  a first  attempt  to  incorporate 
the  frictional  thrust  and  torque  in  the  program;  in  the  event  of  any 
improvement  of  the  estimated  frictional  coefficient,  this  portion  of 
the  program  can  be  easily  modified. 


!■ 

i 


28 


R- 1869 


B.  Thickness  Effects  (Non-lift inj ) 

With  the  propeller  in  translatory  and  rotational  motion,  an  addi- 
tional flow  disturbance  is  generated  about  the  propeller  blade  section 
due  to  the  blade  thickness  distribution;  the  fluid  is  pushed  equally  to 
both  sides,  influencing  the  pressure  distriDution  on  both  suction  and 
pressure  sides  of  the  blade  without  contributing  to  the  net  blade  load- 
ing. As  noted  previously,  this  effect  will  occur  in  steady  flow  condi- 
tions only  since  the  blades  are  considered  rigid.  The  pressure  distrib- 
ution resulting  from  this  "flow  distortion"  thickness  will  be  studied 
in  this  section  by  means  of  the  "thin  body"  approximation. 

The  velocity  potential  due  to  source-sink  distributions  approximating 
the  symmetrical  thickness  distribution  of  an  N-bladed  propeller  is  given  by 


*(x,r,;ri;0)=  - ^ ~L  // 


M , s) 


dS 


n=  1 S 

1 S 

TT  i.  j j 

n=  1 o 


r 1 1 ( . , s ) k 1 + a - o ' 


a - 


pS?sinG  do  op  (1*8) 

b a a ' ' 


where  the  source  strength  is 


5f  (p,s)  of  (C,p,0  ) 

M(p,s)  = 2U  /l  + a2pz — = 2U  — - — — 

d s of 


the  Descartes  distance  R i$ 


[(x-?)2  + r2  + p2  - 2rp  cos  (0  - <p  +9  )] 

o o n 


1/2 


s is  measured  along  the  chord  of  the  blade  section  and  the  otner  symbols 
are  as  defined  in  section  A (see  Eq.  (29). 

From  Bernoulli's  equation,  the  linearized  pressure  P is 


pTU,r,vo)  = -pfU  (£+  a 


3 

3d 


Hence 

P U2  N IT  3f  U,p,0  ) , . 

PT(x,r,<i)o)  = J //  35  “ [(37+a  R]  6b  /r+^P7  sin6ad0adP 

° (49) 
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If  use  is  made  of  the  expansion 


R it 


i m ( u - $ + 0 ) r \ I 

. ° ° n ' (IK)  el(x"Ukdk 


(IK)  = f I ( ! k ; p ) K ( ! kj r)  for  p < r 

ni  J in  rn 

1 I ( | k i r ) K ( | k ! p)  for  r < p 
^ in  1 m 1 


i . im(0-^+6)a>  • / r\  i, 

Ii.+  a _eL.)  >.=  ± y e ° ° n / (k-am)  (lK)m  e‘ U_C)kdk  (50) 


/ o -J  \ I I v 

(37  + a 7T0  R = 7 1 e 

o m=->" 


r = 0 /a  = (op  - 6u  cos0  )/a 
o D Cl 


3fT(C,P,y0) 


3f_ (p  ,6.  ) 


, p . . o u 
0 smO  a 

b ot 


Also  y e'mJn  =(N  f0r  m = ;-N*  1 = ' 

n=l  (0  for  m ^ UN 


On  substituting  (50)  - (52)  in  Eq.  (49),  the  pressure  becomes 
iNpfU2  7T  3f  (p  ,6  ) 

PT(x’r'$o)  = 2tt2  / / /T+a-p^ 

op  a 

i m ( 0 - i>  ) °>  . , ■, . 

• 7 e ° ° / (k-am) (IK)  e‘ (x"C)kdkdS  dp 

m a 

m=-«  -oo 


With  the  trigonometric  transforms t ions  for  0 and  $ 

o o 


i \ - U 

PT  (x-r-o)  = 2~ 1 


> f (p.e  ) 

' ' £ 0__ 

1 : ~ 39 

--■op  a 


. . im(o.  coso  -6, cose  ) 
■) — t -imaa  b a b a 

-p  e e 


...  , -ik(6^cos^  - Gpcos0  )/a 
\ ikAo/a  b a b a , 

l K-aiTj)  ( I K/  e e dxd6  dP  (54) 

m a 


The  0 -integral  is  as  before  (see  Eq.  (33)  of  Section  A) 
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i 

i 


j 

i 

i 

i 

i 


A((k-am)0p/a)  = / 
o 


afT(p,ea) 


38 

a 


i ( (k-am)o£/a)cos6a  dQ 
e a 


(55) 


which  is  evaluated  in  Eq.  (35). 

As  in  the  case  of  "non-plana''"  thickness  effects,  the  thickness 
distribution  can  be  represented  by  Eq.  (31)  and  the  derivative  with  re- 
spect to  0^  by  Eq.  (32).  After  substituting  (55)  and  folding  the  m-series 
to  m = 0 to  +“,  the  pressure  of  Eq.  ($k)  is  brought  to 


iNPfU2  

PT (x, r ,$Q)  = —jp—  I ATa^p7" 

P 

00  i k (Aa-e^cos^  )/a 

• { / A(kQb/a)k(lK)oe  b a dk 

»-oo 

M 00  i (k-am)  (Ao-6  j"cos$  )/a 

+ £ / ( IK)  [A( (k-am)8p/a) (k-am)e  ' 

m=N,2N  m ° 

+ A((k+am)6p/a)  (k+am)e!  (k+am)  ^b"  /a]dk)dp  (56) 

b 


Let  k - am  = u in  the  first  term  and  k + am = u in  the  second  term  of  the 
second  k-integral.  Then  this  integral  can  be  written  (for  p < r)  as 


" /iA(ue“/a)ueiu(“'9bc°%)/a] 

m=N,2N  -»> 

•[I  (ju+am|p)K  (lu+amjr)  + I (ju-am|p)K  (|n-am|r)]du 
m ' m m ' m 1 

= l 2i  /{lmPart[A(u0Ph/a)e!uU°'0b  cos*a)/a]} 

m=N,2N  o ° 

• u tl  (|u+am[p)Km(ju+am| r)  + 1^(|  u-amjp) K^Cj u-am| r) ]du 
and  Eq.  (52)  becomes  finally 

-Np,U2  

pT(*>r’<i’0)  = —^Z — / A +aZ0‘i 

p 

•/0°imPart[A(ueP/a)eiu(Ao'9b  COS*a)/a] 

D 

O 

M 

■u{!  (up)K(ur)+  y [l  (u  +am|  p)  K ( lu+am j r) 
o o 11  m 

m=N ,2N 

+ I (lu-amlp)K  ( | u-am! r) ] }dudp  (57) 

m m 
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for  p < r.  (if  p > r,  p and  r are  interchanged  in  the  modified  Bessel 
functions.  Note  that  a,p,r,u  are  non-d imens iona 1 i zed  by  propeller  radius 

v> 

It  can  be  proved  that  the  integrand  is  zero  when  u = 0.  When  p r , 
the  p-integral  has  a logarithmic  singularity.  (This  is  shown  in  Appendix 
D)  . Although  a logarithmic  singularity  is  integrable,  nevertheless, 
since  the  integration  is  performed  numerically  rather  than  analytically, 
special  precautions  must  be  taken  in  evaluating  the  integrand  in  the  re- 
gion of  the  singularity.  The  procedure  is  described  in  Appendix  E. 


C.  blade  Pressure  Distributions  Due  to  Loading  and  Thickness  Effects 
on  Each  Blade  Face 


The  total  pressure  distribution  on  each  blade  face  is  a super- 
position of  the  pressure  P^  due  to  the  thickness  effect  which  produces 
no  lift,  on  the  pressures  arising  from  the  effects  of  "non-planar"  thick- 
ness, camber  and  angle  of  attack  of  the  blade  and  of  spatial  nonun i form  I ty 
of  the  inflow  field.  The  last  four  components  contribute  to  the  lift  be- 
cause each  produces  a pressure  difference  AP  between  the  back  and  front 
faces  of  the  blade  surface.  On  the  suction  side  (back  face)  the  pressure 
due  to  the  loading  contributors  is  AP/2;  on  the  pressure  side  (front  face) 
the  pressure  is  -AP/2. 


Since  the  blade  is  considered  to  be  rigid,  the  non-lifting  thick- 
ness effect  will  be  present  only  in  the  steady  state.  In  this  case  (q = 0) 
the  blade  pressure  distribution  is  made  up  as  follows: 


a)  On  the  pressure  side 

,AP 


P<°>  . 

p ' 2 w 


(££)  - (■-£) 

V 2 'c  ' 2 ' f 

b)  On  the  suction  side 


<f>npt  + ^ 


(58) 


, (o) 


(-) 

V 2 ;w 


(o) 


+ <f>c  + <f>f  + <f>npt  + P: 


(59) 


where  the  subscripts,  w,  c,  f and  npt  refer  to  wake,  camber,  flow  angle, 
and  "non-planar"  thickness,  respectively. 

For  the  unsteady  flow  case,  q^O 
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a)  r‘«>  . (60) 

b)  p'q)  - (f)‘q>  (61) 

Having  determined  the  pressure  distribution  on  each  side  of  the 
blade  for  all  frequencies  contributing  significantly,  the  instantaneous 
pressure,  which  is  the  sum  of  the  blade  pressures  due  to  all  frequencies, 
can  be  evaluated  as  the  blade  swings  around  its  shaft  (see  the  following 
section  D) . The  instantaneous  pressure  on  the  front  (pressure)  or  back 
(suction)  face  at  a given  radial  position  r is,  when  the  propeller  swings 
around  its  shaft  in  the  clockwise  direction  from  its  upright  position 
(at  12  o'clock) : 


P (O  - Re  l 

' ’ n — 


(r) 

P.s 


- iqO 
e 


= 1 |P^(r)|cos(q0 

q=o 


where 


* 

q 


= phase  angle  = tan 


(P 


(q) 

P,s 


) 


) 


I m 


Re 


(62) 


and  the  subscr  i pts  Re  and  Im  indicate  the  real  and  imaginary  parts  of 

P^!(r).  (9  and  9 are  negative  in  the  clc  i se  direction.) 

P>s  q 

D.  Information  Relative  to  Blade  Strew  .-.na  lysis  an:  Cavitation  Study 

Having  determined  a realistic  blade  hydrodynamic  loading  at  all 
frequencies  of  importance  for  a propeller  operating  in  non  -ur,  i form  in- 
flow, taking  into  account  the  exact  geometr/  (skewed  or  unskewed  blade 
on  a helicoidal  surface)  arid  the  Mutual  interaction  of  the  blades,  both 
the  structural  analysis  of  tne  blade  (displacement  and  state  of  stress) 
and  the  cavitation  study  can  now  be  undertaken. 

A stress  analysis  based  rn  the  finite  element  approach  and  the 
correspond i ng  program  adaptable  to  high-speed  digital  computers,  both 
developed  by  Mechanics  Research,  Inc.,  are  now  available  at  CDC  computer 
centers  world-wide.  Designated  as  STARDVNE,  this  program  has  been  used 
extensively  for  static  and  dynamic  structural  analysis. 


i 
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Experience  at  Davidson  Laboratory^  with  this  program  for  the  struc- 
tural analysis  of  a propeller  blade  has  suggested  the  format  of  the  out- 
put of  the  DL  program  for  blade  loading  distributions  to  suit  the  required 
input  to  the  STARDYNE  program.  The  DL  program  provides  the  mean  and  time- 
dependent  blade  loadings  at  all  significant  shaft  frequencies  at  19  points 
along  the  chords  of  eight  radial  strips.  The  program  can  be  easily 
changed  to  evaluate  the  chordwise  pressures  at  a greater  number  of  points. 
However,  since  the  number  of  radial  strips  governs  the  size  of  the  kernel 
matrix,  a finer  mesh  in  the  radial  direction  will  be  more  easily  achieved 
by  interpolation  of  the  values  between  adjacent  radial  strips. 

The  instantaneous  blade  loading,  which  is  the  sum  of  the  loadings 
due  to  all  frequencies  with  proper  phasing,  is  applied  as  an  input  to  the 
static  STARDYNE  program.  The  angular  interval  between  successive  blade 
positions  can  be  refined,  depending  on  the  type  of  wake.  In  the  case  of 
a container  ship  with  sharp  wake  peaks  at  12  and  6 o'clock,  a finer  spatial 
mesh  is  required  and  a small  angular  interval,  such  as  10  degrees,  is  ad- 
visable during  the  complete  revolution  of  a blade. 

The  blade  pressure- j ump  distributions  due  to  loading  at  each  shaft 
frequency  can  be  extracted  from  the  output  of  the  main  program  and  recorded 
on  a magnetic  tape  to  be  used  in  an  auxiliary  program  to  evaluate  the  in- 
stantaneous pressure  jumps  in  a format  required  by  the  STARDYNE  program. 

Another  feature  of  the  DL  program  provides  information  for  prediction 
of  cavitation  inception.  The  blade  pressures  due  to  loading  and  thickness 
effects  on  suction  (back)  and  pressure  (face)  sides  of  the  blade  which  have 
been  calculated  at  all  significant  shaft  frequencies  will  be  similarly  ex- 
tracted from  the  general  program  and  recorded  on  another  magnetic  tape. 
Another  auxiliary  program  will  use  this  tape  for  the  evaluation  of  the 
instantaneous  pressure  coefficients, 


CP  " 1 


2 PfU2(l  + a2r2) 


(63) 


on  suction  and  pressure  sides  of  the  blade  separately  as  the  blade  swings 
around  its  shaft  (see  Eq.  (62)).  These  values  must  be  scanned  to  determine 
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the  minimum  pressure  coefficient  for  use  in  predicting  the  locations 
for  cavitation  inception. 

Also  extracted  from  the  general  program  and  recorded  on  magnetic 
tape,  for  possible  future  use,  are  the  loading  solutions  L^’n^(r)  for 
all  shaft  frequencies  q and  mode  shapes  n at  the  various  radial  posi- 
tions r. 


to  determine  the  blade  pressure  distributions  due  to  loading  and  thick- 
ness effects  have  been  adapted  to  a high-speed  digital  computer  (CDC  6600 
or  7600).  The  computer  program  provides  the  resulting  hydrodynamic 
forces  and  moments  and  blade  bending  moments  as  well  as  the  blade  pres- 
sure distributions  on  each  blade  face  in  steady  and  unsteady  flow  condi- 
tions. As  noted,  the  unrealistic  leading  edge  singularity  introduced 
in  the  blade  chordwise  pressure  distribution  by  the  cotangent  term  of 
the  Birnbaum  distribution  is  removed  by  the  Van  Dyke-Lighthi 1 1 correc- 
tion method.  In  addition,  the  frictional  contribution  to  steady-state 
thrust  and  torque  has  been  approximated  by  using  the  Prandtl -Schl icht ing 
formula  for  the  friction  coefficient,  without  taking  into  account  the 
viscous  pressure  drag  or  the  drag  associated  with  the  deviation  of  the 
local  flow  direction  from  that  along  the  nose-tail  line. 

A set  of  calculations  has  been  performed  for  the  series  of  3“bladed 
propellers  of  expanded  area  ratio  EAR  = 0.3,  0.6  and  1.2  which  had  been 
studied  in  Reference  3 under  the  assumption  of  zero  thickness,  and  for 
which  experimental  data  were  available  from  controlled  NSRDC  tests1 5 > 1 6 • 1 7 
in  nonuniform  inflow  with  wake  structure  rigidly  specified  (screen  wake). 
The  propellers,  with  destroyer-type  blade  outlines  and  modified  NACA-66 
section  with  a = 0.8  mean  line,  were  one  foot  in  diameter  with  hub  diameter 
0.2  feet  and  pitch  ratio  1.08  at  0.7  radius.  Relevant  geometric  charac- 
teristics of  this  set  of  3-bladed  propellers  are  given  in  the  tables  below. 
The  ratios  of  maximum  camber  to  chord  length  are  given  in  Table  I,  of 
maximum  thickness  to  chord  length  in  Table  II  and  of  leading  edge  radius 
to  chord  in  Table  III. 

TABLE  I 

Ratio  of  Liax i'niun  Camber'  to  Chord,  m /a 


Radius 

tvj 

ci 

11 

Cq 

EAR  = 0.6 

EAR  = 1.2 

.25 

.0370 

.0228 

.0196 

.35 

.0388 

.0231 

.0202 

■ 45 

.0372 

. 022^ 

.0196 

• 55 

.0340 

.0212 

.0185 

.65 

.0312 

.0203 

.0177 

.75 

.0290 

.0198 

.0170 

.85 

.0270 

.0189 

.0160 

.95 

. 0247 

.017*4 

.0147 
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TABLE  II 


Ratio  of  Maximum  Thickness  to  Chord , t /c 

Radius 

EAR  = 0.3 

EAR  =0.6 

EAR  = 1.2 

• 25 

.253 

.090 

.032 

• 35 

.192 

.068 

.024 

• 45 

.146 

.052 

.018 

• 55 

• 113 

.040 

.014 

• 65 

.087 

.031 

.011 

• 75 

.068 

.024 

.0086 

• 85 

.052 

.018 

.0066 

• 95 

.045 

.016 

.0057 

TABLE 

III 

Ratio  of  Leading-edge 

Radius  to  Chord . P 

0 

Radius 

EAR  = 0.3 

EAR  = 0.6 

EAR  = 1.2 

• 25 

.0410 

.00525 

.00066 

• 35 

.0236 

.00290 

.00039 

• 45 

.0136 

.00170 

.00021 

• 55 

.00817 

.00100 

.00013 

• 65 

.00484 

.00060 

.00008 

• 75 

.00296 

.00035 

.00005 

.85 

.00173 

.00025 

.00003 

• 95 

.00133 

.00020 

.00002 

The  coefficients  a^(p)  in  the  thickness  distribution  given  by 
Eq.  (31)  are  presented  in  Table  IV  for  all  three  propellers  and  all 
radial  positions.  The  thickness  distribution  variations  with  EAR  and 
radial  position  are  determined  through  the  t J c , p^  and  0^  factors  of 
Eq.  (31).  (Note  that  the  procedure  for  calculating  an(p)  is  given  on 
page  18.) 


TABLE  IV 


Coefficients  a (p) 
n 

n aR(p) 


1 

2 

3 

4 


0.5^91 

-0.1165 

0 

0 


The  tests  had  been  conducted  in  the  NSRDC  24-inch  water  tunnel 
using  the  closed-jet  test  section  with  a screen  to  produce  the  wake,  a 
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3- cyc)e  screen  for  the  force  along  and  moment  about  the  longitudinal 
axis  and  a 4-cycle  screen  for  the  forces  along  and  moments  about  the 
transverse  and  vertical  axes.  The  tests  were  run  at  a constant  speed 
of  15  rps  with  free  stream  velocity  close  to  12.5  ft/sec.  The  design 
mean  thrust  coefficient  K^.  was  0.150  (practically  open-water  value  in 

the  3~cycle  screen  wake),  the  advance  ratio  varying  slightly  from  propel- 
ler to  propeller17.  Wake  information  is  given  in  Table  V for  3_  and 

4- cycle  screen  wakes.  Note  that  the  Fourier  coefficients  a , b , A 

q q q 

and  B of  Eq.  (23)f  can  be  easily  determined  from  the  information  given 

q 

in  this  table. 

TABLE  V 

Wake  Information  from  Harmonic  Analysis 
of  3-Cycle  and  4-Cycle  Screen  Wakes  (Reference  (16)) 


Radius 

3-Cycle  Wake 

C d> 

3 *3 

4-Cycle  Wake 
C d> 

4 

0.2 

0.089 

18° 

0.050 

0 

0.3 

0. 186 

10° 

0.134 

0 

0.4 

0.220 

6° 

0.170 

0 

0.5 

0.218 

2° 

0.184 

0 

0.6 

0.203 

0 

0.192 

0 

0.7 

0.212 

0 

0.208 

0 

0.8 

0.230 

0 

0.224 

0 

0.9 

0.252 

0 

0.235 

0 

0.95 

0.251 

0 

0.243 

0 

Longitudinal  Veloci 

ty 

vu  ■ 

C sin (qftt  + $ 
q q 

) = Cqsin 

( — qO 

Tangential  Velocity 

V1 

r/u  . 0 

The  computations  were 

run  for 

a series  of 

advance  ratio 

steady-state  case,  design 

J 

and  J = 

: 0.7,  0.6  and 

0.5,  and 

at  1 

in  the  unsteady  cases.  The 

val ues 

of  design  J were  0.841 

, 0.; 

0.844  for  EAR  * 0.3,  0.6  and  1.2,  respectively. 

The  calculated  mean  thrust  and  torque  coefficients  at  design  J 
which  include  the  effects  of  mean  wake,  camber,  flow  angle  and  nonplanar 
thickness  on  the  loading  and,  in  addition,  the  frictional  thrust  and 
torque,  are  compared  with  experimental  data17  in  Figure  4.  Calculated 
results  are  also  shown  without  the  frictional  part. 
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Figure  5 exhibits  the  steady-state  chordwise  pressure  distribution 
at  one  radial  position,  0.65rQ,  on  the  suction  and  pressure  sides  of  the 
blades  for  all  three  propellers,  at  design  J.  These  were  obtained  from 
Eqs.  (58)  and  (59),  for  the  pressure  and  suction  sides  respectively. 

Figures  6 and  7 compare  the  calculated  blade-frequency  thrust  and 
torque  coefficients  and  phases  with  the  experimental  data  of  Reference 
16  using  the  blade-frequency  (third)  harmonics  of  the  screen  wake  given 
in  Table  V.  The  earlier  calculations  of  Refs.  15  and  3 were  based  on 
wake  measurements  reported  by  DTNSRDC  which  were  taken  with  the  screen 
open  sections  at  a different  position.  In  Figure  6,  the  phases  (electrical) 
are  given  in  accordance  with  the  Davidson  Laboratory  definition  on  page 
25  of  this  report.  Figure  7 shows  the  mechanical  phases  according  to 
Reference  16  and  to  Davidson  Laboratory. 

Figures  8 and  9 present  comparisons  of  calculated  and  experimental 
blade-frequency  transverse  forces  and  moments  of  the  three  propellers 
at  design  J.  The  calculations  use  the  harmonics  given  in  Table  V for 
the  A-cycle  screen  wake  in  which  the  tests  were  conducted1®. 

Figures  10,  11  and  12  show  the  blade-frequency  chordwise  pressure 
distributions  at  0.65rQ  for  the  three  propellers  calculated  as  in  Eqs. 

(60)  and  (61).  The  leading-edge  singularity  has  been  removed  by  the 
Van  Dyke-L i ght h i 1 1 correction  method  (see  Eq.  (40)). 

Calculations  have  also  been  conducted  for  off-design  J = 0.5,  0.6 

and  0.7  for  uniform  inflows.  Figures  13,  1**  and  15  compare  the  chord- 

wise  pressure  distributions  at  one  radial  location,  0.65r  for  the  three 

o 

propellers  at  design  J and  J = 0.5  and  0.7.  These  blade  pressure  distrib- 
utions agree  qualitatively  with  experiments  of  Mavlyudov1®  (USSR).  No 
other  experimental  information  is  available  for  comparison. 

Having  evaluated  the  pressure  distributions  in  the  design  and  off- 
design  flow  conditions,  it  is  possible,  in  principle,  to  determine  the 
cavitation  inception  of  an  operating  propeller.  Cavitation  is  assumed 
to  ensue  when  the  local  pressure  falls  to  the  value  of  the  vapor  pressure 
in  the  fluid.  The  speed  V.(r)  of  cavitation  inception  is  known  to  be 
given  by 


1 


I 


I 

I 

I 
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wr=  fen 


where  V(r)  = U /r+T7?2  = resul tant  velocity  at  a blade  section  at 
radial  location  r 


0 = cavitation  number  or  cavitation  index 

= (total  pressure  - vapor  pressure)/ (y  p^V2) 

and  C = minimum  pressure  coefficient,  p /(lp,V2)  (i.e.,  where  C 
p rm  '2^ f p 

is  most  negative) 


Cav i tat  ion 
long  as  a < 
or  a.  = |C 


s 

-C 


assumed  not  to  occur  when  a > -C  , but  will  persist 

, the  critical  cavitation  indexmbeing  when  V./V  = 
^m  1 


as 

.0 


The  critical  cavitation  indices  have  been  evaluated  as  the  absolute 

values  of  the  most  negative  pressure  coefficients  established  by  the 

theoretical  steady-state  blade  pressure  distribution  curves.  They  are 

plotted  in  Figures  16,  17  and  18  for  various  advance  ratios  versus  radial 

position  and  compared  there  with  the  experimental  observations  reported 

in  Reference  17  for  the  same  set  of  3'bladed  propellers  operating  in 

open  water  conditions.  Since  the  experimental  cavitation  indices17  were 

based  on  speed  of  advance  rather  than  resultant  velocity,  the  theoretical 

pressure  coefficients  have  also  been  evaluated  on  that  base  (C  =p  / (4- p rU2) . 

pm  z t 

Although  no  experimental  values  of  a.  for  design  J are  given  in 
Reference  17,  the  theoretical  values  shown  in  Figs.  16-18  for  design  J 
are  justified  by  the  experimental  trends  (i.e.,  extrapolation)  in  that 
vicinity.  At  off-design  J,  the  theoretical  predictions  are  sometimes 
conservative,  sometimes  non-conservative  depending  on  AJ,  EAR  and  radial 
location . 


Since  inception  speed  is  inversely  proportional  to  the  square  root 

of  the  minimum  pressure  coefficient,  0.  = |C  [,  whenever  the  index  of 

P m 

inception  of  cavitation  a.  of  the  experiment  is  greater  than  0.  of  the 
theoretical  calculation,  i.e. 

o . (exp)  > o . (calc. ) 
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then  U.(exp.)  < U.(calc.)-  This  case  indicates  that  the  experiments 

show  that  cavitation  inception  is  at  a lower  speed  than  that  predicted 

by  the  theoretical  calculations  and  thus  the  prediction  is  non-conservative. 

On  the  other  hand,  when 

o . (exp)  < o (calc.  ) 

or  U.(exp)  > U.(calc-),  the  predicted  values  for  cavitation  inception 
are  conservative. 

figure  16,  which  is  a comparison  of  calculated  and  experimentally 
observed  cavitation  index  for  the  3"bladed  propeller  4132  (EAR  = 0-3). 
shows  conservat  /e  predictions  when  J = 0.7  or  design  J,  whereas  the 
predictions  are  non-conservat i ve,  somewhat  at  J = 0.6  and  especially  so 
at  J = 0.5.  except  near  the  hub  and  near  the  tip.  In  Figures  17  and  18, 
for  the  3"bladed  propellers  4118  (EAR  = 0.6)  and  4133  (EAR  = 1.2),  respec- 
tively, all  predictions  are  on  the  conservative  side  except  those  for  J = 0.5. 

Judging  from  this  set  of  calculations,  the  predicted  values  of 
cavitation  inception  are  on  the  conservative  side  compared  with  experi- 
mental observations  as  long  as  AJ  ~ J , - J , < 0.3.  Nevertheless,  the 

d od 

d i sc repanc  ies  which  exist  between  predicted  values  of  the  index  of  cavita- 
tion inception  and  the  experimental  observations  require  investigation. 

These  observed  d i screpanc i es  must  not  be  solely  attributed  to  the 
lack  of  accuracy  of  the  linearized  theory  wnen  modified  by  the  Van  Dyke 
correction.  It  is  well  known  that  visual  determination  of  inception  of 
cavitation  is  dependent  upon  the  subjective  evaluation  of  the  observer 
and  is  highly  dependent  on  the  surface  finish  and  accuracy  of  the  model 
in  the  immediate  vicinity  of  the  leading  edge  of  the  blades.  The  variant 
patterns  of  the  viscous  flow  at  appreciable  angles  of  attack  can  involve 
laminar  separation  with  vortical  flow  standing  off  the  blade  surface 
giving  extremely  low  pressures  in  the  core  of  the  vortex  (as  cited  fre- 
quently by  Eisenberg,  for  example)  and,  hence,  inception  at  o values 
larger  than  the  minimum  pressure  coefficient  provided  by  inviscid  theory 
on  the  surface  of  the  blade.  Also,  from  the  definitive  measurements  of 
Huang  and  Hannan19,  it  can  be  conjectured  that  large  pressure  fluctuations 
can  occur  on  the  blade  surface  at  the  reattachment  zone  abaft  laminar 
separation.  In  addition,  it  is  never  certain  that  observed  inception 
is  attributed  to  truly  vaporous  cavitation  or  to  the  expansion  of  un- 
d i ssol ved  air. 
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The  fact  that  the  experimental  o-values  are  generally  higner  than 

the  cal  cu  1 a ted  I C j at  the  maximum  excurs  ion  in  J (i.e.,  J = 0.5,6J  = 0.33) 

P 

suggests  that  the  real  fluid  effects  alluded  to  above  may  be  responsible. 
However,  at  the  larger  J = 0.7  (smaller  AJ  = 0.13),  the  disagreement  is 
in  the  opposite  direction  except  for  the  first  propeller,  EAR  = 0.30, 
where  the  agreement  is  fine.  Here  the  fluid  speeds  are  of  the  order  of 
twice  as  great  as  at  J = 0.5  and  the  pressure  peaks  on  the  blades  at  re- 
duced angle  of  attack,  while  lower  in  magnitude,  are  indeed  sharper,  i.e., 
the  chordwise  width  of  the  pressure  spike  is  reduced.  One  may  speculate 
that  the  passage  time  of  nuclei  through  the  region  of  low  pressure  is  con- 
siderably less  at  J - 0.7  than  at  J = 0.5  (higher  speed,  shorter  extent) 
and  that  cavitation  does  not  ensue  until  the  a values  are  dropped  below  the 


|C  '• 

Pm  , . 

The  validity  of  the  theoretical  curves  should  be  checked  by  employing 

a rational  method  for  finding  the  effective  angle  of  attack,  the  effective 

camber  and  the  effective  thickness  of  the  two-dimensional  section  (to  be 

substituted  for  each  bl ade sect  ion)  for  which  exact  steady-state  pressure 

distributions  have  been  computed.  This  is  not  a simple  process  as  was 


d i sco ve  red  when  at  temp  ting  to  use  the  families  of  curves  provided  by 
Brockett2®  and  a rule  for  finding  the  effective  angle  of  attack  attributed 
to  W.  Morgan.  The  effective  angle  of  attack  -was  found  to  be  too  impre- 
cise; a very  small  error  produces  a larger  ■„  ror  in  the  C value. 

Clearly,  further  studies  of  the  pressure  distribution  at  the  leading 
edge  of  propellers  are  necessary.  The  results  obtained  thus  far  should  be 
regarded  as  quite  reasonable.  Further  work  may  require  a precise  mechan 
ization  of  the  effective  two-dimensional  solutions  using  certain  inputs 
from  the  three-dimensional  propeller  flow  field,  and  local  sectional 
load  i ngs . 
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CONCLUSION 

In  this  study,  a theoretical  approach  is  evolved  and  a computational 
procedure  adaptable  to  a high-speed  digital  computer  (CDC  6600  and  7600) 
is  developed  for  evaluation  of  the  linearized  pressure  distribution  on 
each  side  of  a marine  propeller  blade,  with  the  objective  to  obtain  suf- 
ficient and  reliable  information  for  cavitation  and  blade  stress  analyses. 

The  essential  information  for  the  blade  stress  analysis  is  the  anti- 
symmetric part  of  the  blade  pressure  distribution,  attributed  to  the  non- 
uniform  inflow  (wake),  blade  camber,  incident  flow  angle  and  nonplanar 
thickness  of  the  blade,  all  of  which  are  associated  with  the  lifting  ac- 
tion of  the  blade  and  contribute  to  the  hydrodynamic  forces  and  moments 
and  blade  bending  moments.  The  symmetric  part  of  the  blade  pressure  dis- 
tribution, attributed  to  the  planar  thickness  and  associated  with  the 
non-lifting  properties  of  the  blade,  contributes  to  the  pressure  dis- 
tribution on  each  side  of  the  blade,  which  is  essential  for  the  predic- 
tion of  cavitation  inception. 

In  addition  to  the  blade  pressure  distributions  at  each  frequency, 
the  program  furnishes  the  propeller-generated  steady-state  and  time- 
dependent  hydrodynamic  forces  and  moments  and  the  blade  bending  moment 
about  the  face-pitch  line  at  the  midpoint  of  any  of  the  radial  strips 
into  which  the  span  is  divided.  The  program  also  provides  the  instan- 
taneous blade  pressure  due  to  loading  alone  and  the  instantaneous  pres- 
sures on  suction  and  pressure  sides,  as  well  as  the  instantaneous  blade 
bending  moments,  as  the  propeller  swings  around  its  shaft.  The  thrust 
and  torque  due  to  friction  is  estimated  by  using  an  approximate  frictional 
drag  coefficient;  in  the  event  of  any  improvement  of  the  estimated  fric- 
tional coefficient,  this  portion  of  the  program  can  be  easily  modified. 

In  applying  the  mode  approach  to  the  solution  of  the  surface  integral 
equation,  the  analysis  and  the  program  are  divided  into  two  main  parts, 
one  dealing  with  the  steady-state  flow  case  (q = 0)  and  the  other  with  the 
unsteady  flow  condition.  The  steady-state  case  is  subdivided  into  design 
and  off-design  conditions.  The  selection  of  the  proper  chordwise  modes 
in  the  steady-state  flow  condition  at  design  advance  ratio  J is  dictated 
by  the  shape  of  the  loading  distribution  in  two-dimensional  flow  on  a 
foil  with  the  same  camber  distribution.  At  off-design  J in  the  steady- 
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state  condition,  there  is  an  angle  of  attack  due  to  the  difference  LJ 
(between  design  and  off-design  advance  ratio)  and  the  additional  load- 
ing due  to  this  angle  of  attack  is  represented  by  the  cotangent  term  of 
the  Birnbaum  modes.  In  the  unsteady  flow  condition,  the  complete  Birnbaum 
modes  are  used. 

The  linearized  unsteady  lifting-surface  theory  requires  the  leading- 
edge  singularity  arising  from  the  cotangent  term  of  the  Birnbaum  modes. 

The  "square-root"  singularity  is  integrable,  but  its  presence  in  the  blade 
pressure  distributions  is  unrealistic  and  has  been  removed  by  employing 
the  Van  Dyke-L ighth  i 1 I correction  factor. 

A set  of  computations  has  been  performed  for  the  series  of  3_bladed 
propellers  for  which  experimental  data  were  available  from  NSRDC  tests  in 
open  water  and  in  the  non-uniform  inflows  due  to  3-cycle  and  4-cycle  screen 
wakes.  The  calculated  results  for  the  hydrodynamic  forces  and  moments, 
steady  and  unsteady,  compare  well  with  the  experimental  values.  There  is 
no  experimental  information  on  blade  pressure  distributions  for  this  set 
of  propellers.  However,  the  blade  pressure  distribution  curves  in  the 
steady-state  case  agree  qualitatively  with  experimental  curves  'l  ®shown  by 
Mavlyudov  (USSR)  for  a different  propeller  model  (NACA-16,  a = 0.8  mean 
line  section,  EAR  = 0.95)  at  0.8  radius.  In  the  absence  of  experimentally 
measured  blade  pressure  distributions  for  the  propellers  treated,  a compar- 
ison is  made  indirectly  through  the  index  of  cavitation  inception  in  uni- 
form inflow.  It  is  seen  that  the  predicted  values  of  a.  = |C  ! are  con- 

i p 

servative  except  at  the  smallest  off-design  J = 0.5  (largest  <^J  = 0-33). 

The  cavitation  index  |CD  j reflects  the  blade  pressure  in  the  neigh- 
borhood  of  the  leading  edge  and  this  is  dependent  on  the  correction  method 
for  removing  the  leading  edge  singularity  of  the  theoretical  distribution. 
However,  the  observed  discrepancies  between  theory  and  experiment  cannot 
be  attributed  solely  to  the  lack  of  accuracy  of  the  linearized  theory  when 
modified  by  the  Van  Dyke  correction.  Experimental  determination  of  the 
inception  of  cavitation  is  dependent  on  the  subjective  evaluation  of  the 
observer.  It  is  also  dependent  on  the  surface  finish  and  accuracy  of  the 
propeller  model  in  the  immediate  vicinity  of  the  leading  edge  of  the  blades 
and  on  the  undissolved  air  content  of  the  fluid  as  well.  There  is  a 
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possibility  that  flow  separation  and  vortex  generation  occur  near  the 
f ieading  edge  and  this  is  not  taken  into  account  by  the  theory.  Further 

studies,  both  theoretical  and  experimental,  are  necessary.  The  results 
I obtained  thus  far  may  be  considered  reasonable. 
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FIG. 2.  DEFINITIONS  OF  ANGULAR  MEASURES 


NOTE.  THE  BLADE  REFERENCE  LINE  IS  THAT  CONNECTING  THE  SHAFT  CENTER  WITH 
THE  MIDPOINT  OF  THE  CHORD  AT  THE  HUB 


STEADY-STATE  CHORDWISE  PRESSURE  DISTRIBUTION  AT  0.65 
RADIUS  ON  DTNSRDC  3-8LA DED  PROPELLERS  AT  DESIGN  J. 


EXPANDED  AREA  RATIO  EXPANDED  AREA  RATIO 
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FOR  DL  PHASE  DEFINITION  SEE  FIG.  7.  (WHERE  <p/3  IS  SHOWN  FOR  EAR-0.6) 
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■REF.  16  DEFINITION  OF  PHASE 


FIG. 7.  MEASURED  AND  CALCULATED  BLADE  FREQUENCY  THRUST  AND  TORQUE 
COEFFICIENTS  FOR  DTNSRDC  PROPELLER  4118,  EAR=0.6, USING  WAKE 
VALUES  OF  TABLE  V FROM  3- CYCLE  SCREEN 


i 


DESIGN  J 


BLADE-FREQUENCY  TRANSVERSE  FORCES  AND  MOMENTS  OF  DTNSRDC  3-BLADED  PROPELLERS 

(IN  4-CYCLE  SCREEN  WAKE) 


DESIGN 


BLADE-FREQUENCY  CHORDWISE  PRESSURE  DISTRIBUTION  ON  EACH  SIDE  OF  PROPELLER  4 1 1 8 
(EAR=0.6  ) AT  0.65  RADIUS 


DESIGN 


BLAOE- FREQUENCY  CHORDWISE  PRESSURE  DISTRIBUTION  ON  EACH  SIDE  OF  PROPELLER  4133 
(EAR*  1.2)  AT  0.65  RADIUS 


= P/l  f u ) ( 1 + Q 


FIG.  13.  STEADY-STATE  CHORDWISE  PRESSURE  DISTRIBUTION  AT 
0.65  RADIUS  ON  PROPELLER  4132  (EAR  = 0.3)  AT  VARIOUS  J. 


FIG.  14.  STEADY-STATE  CHORDWISE  PRESSURE  DISTRIBUTION  AT  0.65 
RADIUS  ON  PROPELLER  4118  (EAR=0.6)  AT  VARIOUS  J. 
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APPENDIX  A 


Evaluation  of  the  and  8 - Integrals  of  the  Integral  Equation  ( 1 9 ) 


0 l(")<y>  -i  ; 4 «e ' >'“»=>  d3 
o 


TT 


0)r,  _ i 


(y)  = ~ , (l-coscp)  e1  C0SV  dy  = J (y)  - iJ.(y) 


TT  J 
O 


(2) 


(y)  = “ J ( 1 +2coscp)  e ' yCOs:?  d9  = J (y)  + i 2 J (y) 


TI  0 
O 


(m>2) 


(y)  = “ j cos  (m  - 1 )p  elyC0S^  dcp  = i™  1 J_  (y) 


TT  J 
O 


m- 1 


where  J (y)  is  the  Bessel  function  of  the  first  kind  of  order  n and 
n ' 

argunent  y 


, ^ t -izcosG 

2)  A n (z)  = j S(n)  e 


s in»  d9 


a)  Birnbaum  distribution 
/ \ ft 

. ( 1 ) , v L ' -3  * 1 zcosS  . . , . , r \ ■ , t \ 

A (z)  * “ j cot  - e s m-  d-:  = J (z)  - i Jj  (z) 

o 

, (n>l)/  \ 1 .■>  • • - -izcosS 

A (z)  = — sin(n-l).'  sintj  e dG 

TT 


/ • »n-2  _ 

L-Ll 1 


J_  (z)  + J_(z) 
n-2  n 


b)  "Roof-ten"  distribution  (a  mean  1 ines) 

/ , \ cos  1 (I  -2a) 

, (1)  , . ' -izcosG  . . , . 

A ( z ) = e sinu  ds 


cos  ' ( I -2a) 


- ' i z J i , 

= e < T + 


(l+COSp)  - I ZCOSd  . , , , 

— =\  e sint;  d 3 

2(1 -a) 


z 2 ( 1 - a)  z; 


| i2az  i2zl  \ 

- Lc  - c I j 


1 


p 
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APPENDIX  3 


TABLE  OF  INTEGRALS  OF  EQUATION  (ZZ) 


= i>  cos9  s inr.d9  = ~ sinX  = 2G(a.)  (=2  when  K = 0) 


. T si/cosu  • ,q  _ 1 c ; 

I)  , e sm-d9  si 

o 


2)  , e 


— I A.  CO  S A* 


s i n2  °dG 


_ iti  , SJJB.  . COS‘  = =4i  F (>. ) (=0  when  X=0) 


- X X X 


3)  ; e 


ri^COS'  s i n 3 Od 0 = f {sinX^3-  + f cos>' ; = 6G(>0  ‘ ^ F(^ 


(=2/3  when  X=0) 


4)  ' e-i>  cos-  sjn40cjo  = ± 


81  (sjnX  J2  + cos^  ,11  . \ y 

A \ A V ,2/  x 2 


~ |4G(x)  + (x-  “ F(>)'  (=0  when  X=0) 


(-D  J, 


5)  ; g-'X.cosS  cos  | d0  = 2Jo(X)  - 4 (4n  + l ) (4n- 1 ) + 4'  n=,  (4n-l)  (4n-3) 


_ B(±x)  (=2  when  X=0) 


(Three  or  four  terms  of  this  series  are  sufficient-) 


i 


{ 


I 


I 
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APPENDIX  C 

RESOLUTION  OF  FORCES  AND  MOVENTS 


With  the  present  coordinate  system  and  sign  convention  (see 
Figures  1 to  3)  the  propeller,  its  N blades  rotating  with  angular 
velocity  -?i  , is  assuned  to  lie  on  a helicoida)  surface  given  by 


F(x,y,z)  = x + ~ tan  — = 0 (C-l) 

3 Z 

where  in  cylindrical  coordinates 

x = cp  /a 

y = -r  s in  6 

z = r cos  & 

9 = cp  - fit  + 8 
o n 


The  unit  normal  to  this  surface  has  components 
F , F , F 


X V z 

n = /8 a^—  5- 

„/F  + F + F 

X y Z 


s a 

1 , z/a  r , ~v/ a r 
/ — zs~ 

+a  r /ar 


so  that 


ar 

“ V = cos  8 
+ a r 


z/r 


n = 


y ~ J\  + aV 


= s i n 3 cos  9 


- y / r 

n = = sin  3 „ j n 0 

Z J)  + a r 

where  3 = tan  ' — , the  hydrodynamic  pitch  angle. 
3 r 


(C-2) 


I 
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The  elemental  forces  are  then 


AFx  = &P  cos  6 AS 


AF  = aF  sin  3 cos  (Qt-C?  -0  ) AS 

y on 

AFz  = - A?  sin  3 s in (Qt-9Q-6n)  AS 


The  elemental  moments  can  be  expressed  as 
i j k 

AQ  = x y z (AP)  (A  S) 

n n n 

x y z 

so  that 


AQ  = A P (yn  - zn  ) AS  = - AP  rsin3  AS 
x z y 


AQy  = - AP(xn_,-znx)  AS 

= AP-r  tanS  sing  s i n (Qt-^0'  6n)  + cosB  cos 


(Qt-tPo^n). 


AQ  = AP(xn  - yn  ) AS 
z y x 

= AP-r  1 cp  tanB  sing  cos(flt-cp  ~e  ) - cosg  s i n (flt-cp  _ -6  ) AS 
o on  o n j 

The  total  force  in  the  x-direction  (thrust)  is 

N iq(Qt-fi)  / n 

F = Re  E e j j AP  q (r.cp  ) cosB(r)  dS 

* n-1  S 

Since  dS  - rdr  dcp  = r%[  sincp  do  dr,  0 s a S ri 

O D Cf  Cf  0( 

L<'C,'>  ^ r ’ = AP^’1  ^r,9o^  r 6b 

N ±iqun  r N f°r  9 = ^N>  £=0,1,2,-.- 

and  Z e = j 

n=l  k 0 otherwise 
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F = Re  •,  N ei;'N_t  ' L ^ ( r ,y  ) cos  3 ( r ) siny  dc  dr 


But  \ L^N^(r,c  ) since  dy  = L^^(r)  the  spanwise  loading 
J a a u 


the  re  fore 


F = Re  Nr  e'/N"t  \ L^'^(r)  cos3(r)  dr' 

x o “ 

The  total  Force  in  the  y-direction  is 
N iq(.Tt-;*  ) / \ 

F = Re  I e " j j £-P  q ( r , c;  ) s i n =(r)  cos  (wt-«pQ- 9p)  dS 

y n=  1 S 

-Re  “ i •L(«)(r.s),inS(r)[.,<^><-'Jt-3.>.,S;“W» 

. 2 i>  c -j  u 

n=]  or 

i (q-l)  (:.t-"„)  -i  l cose,  -i 

+ e e b ' sin  3 d?  dr 

J or  a 

N ±i (q=l ) vp  f N for  q =1  = 2N 
Since  e ‘ l 0 otherwise 


■n  u ia^cos--;  UN+1)  -i  9 cost; 

y ' or 


sin  3 (r)  sine:  dts  d r - 

cc  a 


The  i n tegra  1 s are 


' (T.N+1 ) , b C°S'"a  . 

L (r  ,OrJ  e sin  ^ 


TT  _ . ( ■ N+l ,n)  , v ,-s 
= j _Z  !-  (r)  t(n)e 


;i-3l  cost; 
b a 


s i nt;  dec 
a a 


= _I  L<««?l,S)(r)  A*"1  (ioj) 


where  /3n  (+9^)  is  given  in  Appendix  A. 
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The  re  fore 


Fy  = Re<  — e 


, Nr  i £Nwt  1 

o 


I [' 

o n=  1 


UN-1, n)  (n) 


(r)  A (--/b)  + L 


(iN+!  ,n)  (n) 

(r)  t.  0 


s in  0 (r)  dr 


; 


The  vertical  force  is 


, N i q (fit- 8 ) ^ / -,  . 

F = Re  Z e R I AP'q '(r,t?  ) sing(r)  sin(CTt-0  -0  ) dS - 

2 , J o on 


is  finally 

,-Nr  i JlNQt 
F = Rc  — 2 e 

2 . 2 I 


J 

o n= 


I L 


[' 


(jMI-l.n)  (n) 


sin  g(r)  drj 


n iq(nt-ej  . „ (q) 

Q = Re  J-  I e 


n=  I 


n p ? 

J j 

s 


development  for  Fy 

, the  fore 

(n)  (1N+I ,n) 

(n) 

A (-8J)  - L 

(r)  A (5j 

(C-8) 

')  is 

> 

(r,'vo)  s ing(r)  r dS  - 

and  by  analogy  with  F^  this  becomes 


r 2 ij&Nflt  J UN) 

Q = Re  ■!■  Nr  e f L (r)  s i ng  ( r)  r dr  f 

x o J J 


(C-9) 


The  bending  moment  about  the  y-axis  is 
N iq(nt-0  ) A 

n A 


Q.  = Re(  I e 
^ n=1 


o n 


n J J & p ^ (r,c?Q)  [ rcpQ  tang(r)  sin3(r)  s i n U t-'-?„"9  J 
S 

+ r cosp(r)  cos (Qt-p  -Rn) ] ds[ 


; ch  with  the  trigonometric  transformations  employed  before  becomes 
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- N , , . 2„  / \ r i(q+l)(.it~o  ) mj,  cosy 

„ „'r-  I • • . (q)  , x ■ r sin  - (r)  r n'  b cc 

Q - Re  Z t i i L (r,y  ) m,  cosy  e e 

y n=l  2 J J v b cosp(r)  b 


i (q- I ) (fit-0  ) - i o!"  cosy  i (q-i  i ) (.it-5  ) i gT  cosy 

n b c/ 1 ri  d Of 

- e e J + cos  3 ( r) 1 c e 

i(q-l)(at-0  ) - I pf  cosy 

+ e e r sin;  d;  dr 

J j j 

N iANZt  r n (;.N- 1 ) i choosy 

= Re  - e 'j*  1 | ip"  tanj3(r)  sing(r)  L (r,y  )e  ' 

r O 


(XN+l) 

- L (r,y^)  e 


- i 0,  cosy  - 
b To  | 


• cosy  siny  dy 

u'  Ol  0 


-t  (xN-l)  iPy  cosy  (XN+1)  -i9*"  cosy  - 

+ j'  cos3(r)j_L  (r,ya)e  b a+L  (r.y^e  b ®! 


siny  dy  • ' r dr 


The  f i rs t y - in tegral  is 

a 

~ UN+1)  = ioj*  cosy 

,•>  , i \ b a . 

i L (r,y  ) e cosy  s my  dy 

*)  a a a ot 


n (TN+1 ,n) 

i Z L 


*j  i 

o n=l 


=i  -S  cosy_ 

(r)  o(n)  e " cosy  siny  dy 

a a a 


(/,N+I  ,n) 


= Z L 


(r)  A®  (i.-br) 


(C-10) 


whe  re 


( ) is  as  defined  in  Appendix  A-  The  second  y -integral 


is  given  by  ( C -6) - Finally, 


Nr  i AN.,  t 1 

Q.  = Re  — ~ e ‘ i \ s 

y .2  b 

o 


i n 3(r)  tang(r)  _Z  j L 

l ■— 


UN-l,n)  /-\ 


(TN+1  , n)  /-x  r (AN-1  .n)  (n) 

- L (r)  ^(0b)J  + cos^ ( r ) _Z  |^L  (r)A  (-6b) 

n=l  L 


(AN+1  ,n)  (n) 

+ L (r) 


(C-ll) 
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The  bending  moment  about  the  z-axis  is 
N iq (Qt-Q  ) „ / x f 

Qz  = Re  E e P J | APVq  (r,c?Q)[_rtpo  tanp(r)  sin^(r)  cos  (r<t-<pQ- «n) 


_ 

- r cosg(r)  s in  (fit-cp  - 6 ) dS 


It  can  be  shown  that 

-Nra  i£NAt  1 


(IN- I , n)  (n) 


(r)  /.,  (-0  ) 


(XN+1 , n) 

O 

(n)  -| 

+ L 

(r) 

N (E'b>] 

+ cosr(r) 

1 

(-^N+l  , n) 

(n) 

- L 

(r) 

K>] 

p-  • r d r - 

(iN-1 , n)  (n) 

(r)  A (-ebr) 


- L (r)  A (e£)  j • • r dr  - ( C - 1 2 ) 

In  the  text,  and  in  the  program  as  well,  the  hydrodynamic  pitch 
angle  3(r)  of  the  assumed  helicoidal  surface  is  replaced  by  the  geometric 
pitch  angle  ~ (r)  of  the  actual  propeller- 
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APPENDIX  D 


THE  SINGULARITIES  OF  THE  INTEGRANDS  OF  THE 
BLADE  PRESSURE  DISTRIBUTION 


In  Equation  (49)  for  the  pressure  due  to  thickness  (non- 1 i f t i ng) 
it  is  seen  that  there  is  a singularity  when  p = r,  x = §,  6q  = c?q  and 
9^  = 0 (i-e.,  n = l).  The  singular  part  of  the  pressure  can  be  expressed 
as  (see  Eq  • (53)  with  the  substitutions  0o  = a§  and  c?  = ax): 


pfUa 


df  (p,9a) 
39 


J 1 +a2  p2  J 


J 

- CO 


k(IK)o  eik<x-5)  dk 


♦ Z j’  <IK)m  [(k-«n)ei<k'm)(,,‘5)+  (k«®,)ei(k+m)<*-5)]  dk^dp 
m=  I 

(D-l) 

For  arbitrary  thickness  (see  Eq.  (32)). 

of(p,er/)  y 4 

— rr  « Co(p)  cosy  + 4 Cn(p)  sin  n~ 

a n=l 

where  the  coefficients  are  obtained  as  shown  i n sec t- A,  3 , c • Then  the  trigo- 
nometric transformation 


§=  cos9a)/a 


yields 


cos  9 = 

a 


- a= 


cos 


f-V- 


9°  + o°-  a? 

b 


2 ar- 


sine d9  = d| 

or  a ,,  p 5 


I 


I 

I 
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sin  2 ' d"  = 2sinj  coss  d- 
i .i  y y y 


2a  (c  -a*) d: 


sin  3l;  d r* 
a or 


■ t +4  ^ 1 

(0LP) 


Cl  - 3 P e \ -| 

? - ai  i . c - a?  a 

. . 4 , 


sin  4=  = [~  8 -4'.^ 

o-  L \ „ p / v c P 

"b  6b 


and  thus  the  slope  can  be  expressed  in  the  following  form: 

(p.e^dA^  ~ |[dQ(p)  + ^dj  (p)  - d2(p)§  + d3(p)§  + d^(p)?  +d ^ ( . )i 
Of 

_ P(?.o) 


Equation  (D-l)  can  be  written  as 

-4-  J J F(«,c)  ^ _i  k(IK)o  e,k<X"-4k 


. ' (IK)  !‘(k-an)4(k'a")(x'd  + (k*8.  iJWMl  dk-  d?ip 

11=1  -io 


For  finite  tii  i Allans  i on  of  1/R  in  the  above  has  no  singularity. 

The  singular  behavior  is  present  only  in  the  infinite  m-series  (see  Ref- 
erence 10).  When  n * M large,  the  generalized  mean  value  theorem  can  be 


used : 


d d 

' f (k)  p(k)  dk  ■-  f (A)  . p (k)  dk,  c 4 A s-  d 


where  f(k)  = I 1^(1  kip)  K^(fklr)  for  p < r 


I ( I k I r)  K ( I k I p)  for  r < o 
m m 


and  f (A)  = 1 ( I A I p)  K (lAIr),  etc-  w i th  A < < m (orde  r)  • 

By  using  N i cho I son  1 s ^ ^ approx imat i on  for  the  product  of  the  modified 
Bessel  functions  when  A < < m 


f(A)  rsT-  Zm  where  Z = ( p/r  for  p < r 

zm  ) 

( r/p  for  r < j 

Then  for  large  m the  integral  can  be  written  as 

ik(x-g) 


l=jjF(5.p)  E - Z -'cos  am  (x-§)  j ke 

5 p m=M  m -so 


+ i am  sin  am  (x-g)  , e 


i k(x-§) 


From  Jolley's  collection  of  series  summations 

iaZ  sin  a (x-r) 


dk 

dk  • d|  d: 

(12) 


(0-3) 


(D-4) 


ia  Z Z sin  ma (x-§)  = 

m=l  1-2Z  cos  a(x-g)+Z 


(see  Jolley  499) 


Also  Z — Zm  cos  ma(x-?)  = — Z Zn  cos  ma(x-c) 


dZ  m l m 

m=  I 


m=  1 


cos  a (x-Q  - 7 
1-2Z  cos  a(x-S)+Z 


~ (see  Jol ley  500) 


The  re  fo  re 


Z — Z™  cos  ma(x-§)  = 


’ cos  a 


(x-S)  - Z 


I-2Z  cos  a(x-§)+Z 


dz 


J 2 _ 

- - log  [1  - 2Z  cos  a (x-g)+Z  J 


The  m-series  of  (D-4)  then  is  equivalent  to 


r I r , . 2,  ik(x-C) 

S = J - — log  [ 1 -2zcosa (x-^) +Z  ] i ke  dk 


iaZsin  a(x-^)  ? 

2 J 

1-2  Zcos  a(x-?)+Z 


ik(x-£)  "■ 

e * dki- 


E — Zm  i cos  am(x-£)  j ke'^X  ^ 


+ i am  sin  am(x-?)  J 


i k(x-§)  dk' 


(D-5) 


where  the  finite  m-series  can  be  ignored  since  it  is  certainly  not  singular- 
The  k-integrals  are  evaluated  as 


J eik(x*?)dk  = 2n  6(x-g) 

“CO 


J k eik(x"5)  dk  = - i 2tt  6 ' (x-g) 
-00 


(D-6) 


where  s(x-?;)  is  the  Dirac  delta  function  and  6'(x-g)  the  derivative  of  this 
function  with  respect  to  (x-g)- 

With  the  substitution  of  (D-5)  and  (D-6)  and  letting  xq  = x - g, 
Equation  (D-^)  becomes 

1 = i2n  J J f(x"-><0«p)  2 '°9  (1_2Z  cos  a*0+z  ^ 6 ' (*Q) 


*0  P 


aZ  sin  ax 

9 5-  6(x  ) r dx  d 

1-2Z  cos  ax+Z  ° ' ° P 


and  integrating  over  x results  in 
, ^F(x'x0)p) 

* ~ +'n  J dx  x =0 

p 00 


log  (1-2Z+Z  ) dp 


(D-8) 
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Wi th  z = p/r  or  r/p , 

1 og ( I -21+Z  ) = log(l-Z)2=  2 1 og ( 1 -2) 

= 2 log  or  2 log  C2-^-) 

hence  the  p-integral  has  a logarithmic  singularity  when  p - r and  this 
has  a finite  contribution  and  is  integrable- 
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APPENDIX  E 

EVALUATION  OF  THE  p- INTEGRATION  IN 
THE  REGION  OF  THE  SINGULARITY  (SECTION  3) 

Let  K(p)  represent  the  integrand  of  the  p- i n teg ra t ion . It  has  been 

seen  that  when  p -•  r,  K(p)  varies  as  in(p-r).  A logarithmic  singularity 
is  integrable,  but  since  the  integration  is  performed  numerically,  special 
precautions  must  be  taken.  In  the  region  of  the  singularity  p = r,  the 
integral  is  put  in  the  form 

r+6 

' M(P)  , .X 


where  M(p)  = (o-r)R(p)  so  that  M = 0 when  p = r,  and  3 = Ar/2. 


The  function  M(o)  can  be  expanded  about  the  singularity  p = r 
by  the  Lagrange  formula 


M(p)  = 


i = 0 , i n 


whe  re 


nn (p)  = (p-oo)  (p-p, ) — (p-Pn) 


n 1 (p  • ) = “7”^  (p)  eval  uated  at  p = o . 
n i dp  n i 


and  M.+|  = H(p.)  (see  Scarborough^  and  Watkins  et  al  1 '). 

In  the  strip  from  r - 3 to  r + (3  (with  n = 4 for  the  5"Point 
formula),  po=r-j3=r_  26,  pj  = r - 5,  etc-  where  6 = 3/2-  Then 

n’(Pj)  = (-1)4"1  ^ il  (4-i)l  (E-3) 


M(p)  = I 

FT  i»  1 1(4-1)  I 


(o-r)  (p-r-c) (p-r-2*) 
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E-2 


where  = 0 since  = r- 
The  integral  is 
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